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warm inflation are discussed. Quantum fleld theory real time calculations at finite temperature 
are then presented and the derivation of dissipation and stochastic fluctuations are shown from a 
general perspective. Specific results are given of dissipation coefficients for a variety of quantum 
—.^ ■ field theory interaction structures relevant to warm inflation, in a form that can readily be used by 

^sL I model builders. Different particle physics models realising warm inflation are presented along with 

, their observational predictions. 

^ ■ PACS numbers: 98.80.Cq, ll.lO.Wx, 12.60.Jv 

a 
^_ 

keywords: cosmology, inflation, dissipation, particle production, supersymmetry 

00 ; 

I In Press Reports on Progress in Physics, 2009 



(N 
> 

in 



X 



I. INTRODUCTION 



It has been fourteen years since warm inflation was introduced and with it the first and stifl only alternative 
(-H ' dynamical realisation of inflation to the standard scenario. The standard picture of inflation introduced in 1981 
relied on a scalar field, called the infiaton, which during inflation was assumed to have no interaction with any other 
fields. During inflation, this field rolls down its potential and due to it being coupled to the background metric, a 
damping-like term is present which slows down its motion. As this infiaton field was assumed to not interact with 
other fields, there was no possibility for radiation to be produced during inflation, thus leading to a thermodynamically 
' supercooled phase of the Universe during inflation. Getting out of this inflation phase and putting the Universe into a 
00 ' radiation dominated phase was a key issue, termed the "graceful exit" problem [1, 2, 3]. The first successful solution 
of the graceful exit problem and so the first successful cold infiation model, was new infiation The solution was 
QQ , to picture particle production as a distinct separate stage after infiationary expansion in a period called reheating. 

In this phase, couplings to other fields were assumed to be present and the infiaton would find itself in a very steep 
00 , potential well in which it would oscillate. These oscillations would lead to a radiative production of particles. There 
' have been many variants of the original cold inflation picture, first introduced in the context of the new inflation 
model and shortly afterwards in the chaotic inflation model Q, with many other models that followed. 

The warm inflation picture differs from the cold inflation picture in that there is no separate reheating phase in the 
former, and rather radiation production occurs concurrently with inflationary expansion. The constraints by General 
Relativity for realising an inflationary phase simply require that the vacuum energy density dominates and so this 
does not rule out the possibility that there is still a substantial radiation energy density present during inflation. Thus 
on basic principles, the most general picture of inflation accommodates a radiation energy density component. The 
presence of radiation during inflation implies the inflationary phase could smoothly end into a radiation dominated 
phase without a distinctively separate reheating phase, by the simple process of the vacuum energy falling faster than 
the radiation energy, so that at some point a smooth crossover occurs. This is the warm inflation solution to the 
graceful exit problem. 

Dynamically warm inflation is realised if the inflaton were interacting with other fields during the infiation phase. In 
fact in any realistic model of infiation, the infiaton must be coupled to other fields, since eventually the inflaton must 
release its vacuum energy to other fields thereby creating particles which form the subsequent radiation dominated 
era in the Universe. Thus the idea that these couplings to other fields somehow are inactive during inflation, as 
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pictured in the cold inflation picture, is something that does require verifying by detailed calculation. When such 
calculations are done, the result is that there are regimes in which particle production during inflation occurs. This 
review will present the calculations which demonstrate particle production during inflation, thereby leading to a warm 
inflationary expansion. 

The idea of particle production concurrent with inflationary expansion was first suggested in the pre-inflation 
inflation paper by L.Z. Fang in 1980 His paper proposed using a scalar field with the origin of infiationary 
expansion due to a claimed anomalous dissipation term that would be generated based on Landau theory if the field 
was undergoing a second order phase transition. This was dynamically very different from the scalar field inflation that 
eventually became successful, and the source of dissipation was also different from that in warm inflation. However 
this model captured the basic idea of concurrent particle production and inflation. Then in the mid-80s two papers 
proposed adding a local T(f> type dissipation term into the evolution equation of the inflaton field. Moss [6| and 
then Yokoyama and Maeda 01 ■ In both cases the dissipative term generated a source of radiation production during 
inflation. The idea of a dissipative term was re-discovered independently by Berera and Fang [0| almost a decade later. 
They went further by proposing that the consistent dynamics of the inflaton fleld was a Langevin equation, in which 
a fluctuation-dissipation theorem would uniquely specify the fluctuations of the inflaton field. That paper by Berera 
and Fang provided the foundations for the theory of fluctuations in warm inflation and the Langevin equation has 
since been the fundamental equation governing inflaton dynamics. Following that work, in Berera proposed that a 
separate reheating phase, as standard in all inflation models up to then, could be eliminated altogether. This paper 
proposed a new picture of inflation, which it termed warm inflation, in which the process of inflationary expansion with 
concurrent radiation production could terminate simply by the radiation energy over-taking the vacuum energy, thus 
going from an inflationary to a radiation dominated era. This work presented an alternative solution to the graceful 
exit problem to the one given by the standard inflation scenario. This warm inflation picture was verified explicitly 
in [lO|, where the Friedmann equations for a Universe consisting of vacuum and radiation energy were studied and 
gave many exact warm inflation solutions to the graceful exit problem. Finally in [flj the calculation of fluctuations 
were done by Berera for the inflaton evolving by a Langevin equation in a thermal inflationary Universe. 

Alongside the development of the basic scenario, the first principles quantum field theory dynamics of warm in- 
flation was developed. This started in with a quantum mechanical model which demonstrated the origin of the 
fluctuation-dissipation relation in warm inflation. The key step in deriving warm inflation from quantum field theory 
is in realising an overdamped regime for the evolution of the background inflaton field. The initial attempt at this was 
done by Berera, Gleiser and Ramos in (Tsj . In this work, it was proposed that the overdamped evolution should occur 
under adiabatic conditions in which the microscopic dynamical processes operated much faster than all macroscopic 
evolution, in particular the scalar field motion and Hubble expansion. Based on this criteria, a set of consistency 
conditions were formulated in p^ . which would be required for a self-consistent solution. However in [l3| no explicit 
warm inflation solutions were found. This point was further highlighted by Yokoyama and Linde in |14| . in which 
several models were studied from which the conclusions of [l3| were verified. The problem in these early works was 
that dissipation effects were being looked at in a high temperature regime and it proved too difficult to keep finite 
temperature effective potential corrections small, so that the inflaton potential remained relatively flat, and at the 
same time obtain a large dissipative coefficient. One type of model was shown that could realise such requirements [Tsj 
and this was the first quantum field theory model of warm infiation, although it was not a very compelling model. 
Subsequently Berera and Ramos in [16] suggested a solution for getting around the mutual constraints of obtaining 
a large dissipative coefficient and yet small effective potential corrections. The main observation was that supersym- 
metry can cancel local quantum corrections, such as zero temperature corrections to the effective potential, whereas 
temporally non-local quantum effects, such as those that underly the dissipative effects, will not be cancelled. This 
led to |16j proposing a two-stage interaction configuration, in which the inflaton was coupled to heavy "catalyst" 
fields with masses larger than the temperature of the Universe and these fields in turn were coupled to light fields. 
The evolution of the inflaton would induce light particle production via the heavy catalyst fields. Since these heavy 
catalyst fields were basically in their ground state, the quantum corrections associated with them could be cancelled in 
supersymmetric models. The calculation of the low temperature dissipative coefficients for this two-stage mechanism 
were first done by Moss and Xiong [17.] . 

There is an earlier review which covered the basics of the warm inflation scenario (isj . In this review full details will 
be developed of the quantum fleld theory dynamics of warm inflation. This will first start in Sec. |TT]with a summary 
of the warm inflation scenario, including a comparison of it to cold inflation. In Sec. IIIII a basic introduction to 
thermal fleld theory is given including the real time formalism for interacting fleld theories. In Sec. IIVI the effective 
evolution equation of the inflaton fleld is derived, in which all fields it interacts with are integrated out, leading to 
a Langevin type nonconservative equation which contains a dissipative term and a noise force term. The detailed 
properties of these dissipation and fluctuation terms is then studied in Sec. |Vl In addition the physical picture of the 
dissipation effects in warm inflation are discussed. In Sec. IVII the calculations are extended to curved spacetime. In 
Sec. IVIII various particle physics models of warm inflation are presented. Finally Sec. IVIIII presents some concluding 
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remarks and future work being done on the subject. Our conventions are as follows. We use spacetime metrics with 

(H ) signature, and we use natural units for the Planck's constant, Boltzmann's constant and the velocity of 

light fi, = fc = c = 1. 



II. THE WARM INFLATION PICTURE AND ITS MOTIVATIONS 



A. Background equations 

The foundation of cosmology is the cosmological principle which states no observer occupies a preferred position 
in the Universe. This principle implies the Universe must be homogeneous (looks the same from every point) and 
isotropic (looks the same in all directions). To quantify the cosmological principle, in regards the geometry, the metric 
of a space-time which is spatially homogeneous and isotropic is given by the Friedmann-Robertson- Walker (FRW) 
metric, 



1 — kr^ 



' sin^ 9d(p'^ 



(2.1) 



where (r, 0, ip) are spherical-polar coordinates parameterizing the spatial dimensions, t is cosmic time, and fc = 1, — 1, 
describe spaces of constant positive, negative and zero spatial curvature. The most important quantity here is a(t) 
which is the cosmic scale factor and describes the expansion of the Universe. The Hubble parameter, defined by 



H=-, (2.2) 
a 

quantifies how fast the Universe expands. 

The evolution of the scale factor is related to the pressure p and density p by the scale factor acceleration equation, 



d AttG , „ , ,„ 

- = ^(p + 3p), (2.3) 

a 3 

where G is Newton's gravitational constant. Energy conservation is expressed by 

p=-3Hip + p). (2.4) 

Inflation by definition is a phase when the scale factor is growing at an accelerated rate, a > 0, which based on Eq. 
(j2.3p requires p < — p/3, thus for a substance with negative pressure. The spatial part of the metric evolves rapidly 
towards the flat metric fc = 0, and 



= 8nGp . (2.5) 

The most common example of inflationary expansion occurs when the dominant form of matter has equation of state 
Pv — ~Pv ~ constant, which is called vacuum energy. Such an equation of state plugged into Eq (|2.5p leads to an 
exponential scale factor behavior a{t) ~ exp(77t), with H constant. 



B. Inflaton dynamics 



The key problem of inflationary cosmology has been trying to realise inflation from a realistic particle physics 
motivated model. The general observation which has driven this idea is that a scalar field has the necessary equation 
of state needed for inflation. In particular the energy and pressure density of a scalar field are, respectively, given by 



P 
P 



V(0). 
V(0) 



2a2 
6a2 ■ 



(2.6) 
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FIG. 1: Comparison of the cold and warm inflationary pictures [31 . Top graphs show the scalar field evolution and the bottom 
graphs show the vacuum and radiation energy density evolution. 

with the key point being that the potential energy of this field has precisely the equation of state conducive to 
inflation. Thus, the basic idea has always been to somehow get the potential energy of a scalar field to dominate the 
energy density in the Universe, and thereby drive inflation. And then, once enough inflation has occurred, to convert 
the potential energy into radiation and enter into a radiation dominated expansion phase. There are two underlying 
dynamical realisations of inflation, into which all models fall, cold and warm inflation. Both dynamical pictures are 
summarized in Fig. [T]and this subsection will review both pictures. 



1. Cold 


inflation 

















FIG. 2: A quadratic inflationary potential with the inflaton initially starting at some large amplitude. 

This is the standard picture of inflationary dynamics. In the cold inflationary picture, one or more scalar inflaton 
fields are assumed to decouple from everything apart from gravity. The energy density is dominated by the scalar 
field potential. For example, in Fig. [5] it is shown a m^^^ potential, in which when the inflaton amplitude is displaced 
to some > 0, inflation can occur. The evolution of the scalar field in the FRW universe is described by the General 
Relativistic version of the Klein-Gordon equation. 



5 



^i + 3i/0-^V2</, + y'(</,) = O. (2.7) 

In this equation the Hubble damping term, 3H(f>, formally acts like a friction term that damps the inflaton evolution. 
However this 3H(j3 term does not lead to dissipative energy production, since its origin is from the coupling of the 
scalar field with the background FRW metric. The inflaton plays the role of driving inflation as well as providing the 
seeds for density fluctuations. This Subsection only focuses on the first of these requirements. 

In order for inflation to occur, the inflaton must be potential energy dominated, which means the potential energy 
V{(l)) must be larger than the gradient energy (V0)^/2 and the kinetic energy 0^/2, 



y(^) » (V^) 2/2, </'V2. (2.8) 

Moreover, in order to obtain enough inflation, these conditions must persist for some span of time. The usual way to 

achieve this is for the inflaton to start out almost homogeneous and at rest in some small patch of space and then to 
have the inflaton evolution equation overdamped, with approximate form 

3i7(/) + y' ((/))« , (2.9) 

« 8TTGV{(j)) . (2.10) 

The consistency of this approximation is governed by conditions on a set of two slow-roll parameters, 



lip 



«1, 



r, = « 1, (2.11) 

where mj,^ = SttG, so mp — 2.4 x lO^'* GeV. If these conditions hold in a region of space, then inflation can happen. 

For the inflaton dynamics described by Eq. (|2.7p . it is instructive to see how an inflationary scale factor growth 
occurs. It follows from the approximate Eqs. (|2.9p and (|2.10[) that H — —e and therefore, for e <C 1 the expansion is 
roughly constant. From Eq. (|2.10p we see that a{t) w a(0) exp{Ht), where w 8ttGV/3, as we stated earlier. 

The evolution of any radiation contribution to the energy density in the universe also can be easily studied in this 
example. The energy conservation equation becomes pr = — 4iJpr, which has the exponentially decaying solution 
Pr ~ pr(0) exp(— 4_fft). In other words, whatever the initial radiation energy density in the universe at the onset of 
cold inflation, this rapidly decays away, thus supercooling the universe. As shown in Fig. [1] during cold inflation, 
the vacuum energy density is large and almost constant, whereas the radiation energy density is negligible, hence the 
name cold inflation. 

Once a region of adequately large potential energy materializes, the physics of the subsequent evolution is quite 
straightforward. The gravitational repulsion caused by the negative pressure drives that region into a period of 
accelerated expansion. One expresses the amount of inflation as the ratio of the scale factor at the end of inflation 
qei to that at the beginning ubi, and it is usually stated in terms of the number of e- folds Ng, 



N, = \n—. (2.12) 

asi 

In order to inflate the horizon size to a scale covering the observable universe, it is necessary to have around 60 e-folds 
of inflation. 

Eventually, inflation must end and radiation must be introduced into a very cold universe so as to put it back into 
a radiation dominated Hot Big Bang regime, which is the so calledgraceful exit problem [ll,[3, S]- In the cold inflation 
picture, the process that performs this task is called reheating [Sj. This is usually envisioned as occurring shortly 
after the slow-roll approximation has broken down, and is often associated with oscillations of the inflaton field about 
the minimum of its potential. An example of this type of potential will be studied in more detail in Subsec. Ill Dl If 
the inflaton is interacting with other matter fields, the oscillations of the inflaton will lead to particle production so 
that, as shown in Fig. [H the radiation energy density begins to increase. 



2. Warm inflation 



The other dynamical reahsation of inflation is warm inflation. This picture is similar to cold inflation to the extent 
that the scalar inflaton field must be potential energy dominated to realise inflation. The difference is, in this picture 
the inflaton is not assumed to be an isolated, non-interacting fleld during the inflation period. So, rather than the 
Universe supercooling during inflation, instead the Universe maintains a small amount of radiation during inflation, 
enough to noticeably alter inflaton dynamics. In particular, the dividing point between warm and cold inflation is 

roughly at pV'^ ~ H , where pr is the radiation energy density present during inflation and H is the Hubble parameter. 
Thus pV^ > H is the warm inflation regime and ~ is the cold inflation regime. This criterion is independent 
of thermalisation, but if such were to occur, one sees this criteria basically amounts to the warm inflation regime 
corresponding to when T > H . This condition is easy to understand since the typical inflaton mass during inflation 
is ~ H and so when T > H , thermal fluctuations of the inflaton fleld will become important. Subsequent to 
the introduction of warm inflation, other scenarios have been suggested which utilize some of its concepts such as 
thermal fluctuations during inflation and graceful exit via radiation energy density that smoothly becomes 

dominant [23|; however the dynamics in these scenarios differs from that of scalar field stochastic evolution. 

The interaction of the inflaton with other flelds implies its effective evolution equation in general will have terms 
representing dissipation of energy out of the inflaton system and into other particles. Berera and Fang ^] initially 
suggested that for a consistent description of an inflaton field that dissipates energy, the inflaton evolution equation 
should be of the Langevin form, in which there is a fluctuation-dissipation relation which uniquely relates the fleld 
fluctuations and energy dissipation. This has formed the basis of all fluctuation calculations in warm inflation. The 
simplest such equation would be one in which the dissipation is temporally local. 



4>+[3H + T]^~-j^V^^ + V'i<f>)=C (2.13) 

In this equation, T(f> is a dissipative term and ^ is a fluctuating force. Both are effective terms arising due to the 
interaction of the infiaton with other fields. In general these two terms will be related through a fluctuation-dissipation 
theorem, which would depend on the statistical state of the system and the microscopic dynamics. Equations like 
(|2.13p are the main subject of this review, and will be further described in the next subsection and in Sees. IllUfVl 

In order for warm inflation to occur the potential energy py must be larger than both the radiation energy density 
Pr and the inflaton's kinetic energy. A major difference from cold inflation is in the evolution of the energy densities, 
as can be compared in Fig. [1] In warm inflation the radiation energy does not vanish because vacuum energy is 
continuously being dissipated at the rate p^ = — T(/>^. The energy conservation equation (|2.4p for this system of 
vacuum and radiation becomes 



Pr = -4Hpr + T</)2. (2.14) 

In this equation the second term on the right-hand-side acts like a source term which is feeding in radiation energy, 
whereas the first term is a sink term that is depleting it away. When H, T and are slowly varying, which is a good 
approximation during infiation, there will be some nonzero steady state point for pr- Thus at large times, compared 
to the local Hubble time, the radiation in the universe becomes independent of initial conditions and depends only 
on the rate at which the source is producing radiation. 

A series of slow-roll conditions must be satisfied for a prolonged period of infiation to take place. The slow-roll 
parameters for warm inflation are 

The slow-roll conditions for warm inflation can be expressed as 



TV 
TV 



(2.15) 



e<l + Q, TjKl + Q, f3<l + Q, (2.16) 
where the parameter Q is defined by 
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The inflationary solution to the system of equations can be shown to be an attractor when the slow-roll conditions 
are satisfied ^] . These conditions can be weaker than the corresponding slow-roll conditions for cold inflation if Q is 
large. Inflation ends when the vacuum energy ceases to dominate, which is typically when p„ = pr and e = 1 + Q. An 
example is shown in Fig. [TJ Exact solutions for background warm inflationary cosmologies of radiation and vacuum 
energy densities were computed in |,10.] . 

Additional slow-roll conditions must be imposed if the dissipation coefficient or the potential depend on the radiation 
density. For example, in the case of thermal radiation there are quantum thermal corrections to the inflaton potential. 
An additional slow-roll parameter describes this effect. 



(2.18) 

The slow-roll condition for S is stronger than the slow-roll conditions on the other parameters [2^, 

S<1. (2.19) 

This condition is crucial to the realisation of warm inflationary models. Basically, this condition states that viable 
warm inflationary models use some mechanism for suppressing thermal corrections to the inflaton potential. 

It is worth stressing the fact that the presence of radiation during the inflationary epoch is perfectly consistent with 
the equations of motion. All that is required is that the vacuum energy density py be larger than the radiation energy 
density p,-. In cold inflation, both the radiation energy density and the friction term are negligible, i.e., p^^ ^ H 
and T ^ H. The main difference between warm inflation and cold inflation is the reversal of this first condition, i.e. 
pl^'^ > H. There are two regimes that can then be identified, strong and weak dissipative warm inflation. Strong 
dissipative warm inflation is the case T > 3H, and weak dissipative warm inflation is T < SH. The terminology here 
is almost self-explanatory, in the strong dissipative regime, the dissipative coefficient T controls the damped evolution 
of the inflaton field and in the weak dissipative regime, the Hubble damping is still the dominant term. 

Even though the presence of radiation need not hinder inflationary growth, it can still influence inflaton dynamics. 
Consider, for example, inflation occurring at the Grand Unified Theory scale, which means 

this case the Hubble parameter turns out to be ^ V^^^/mp ^ 10^° GeV. For cold inflation and weak dissipative 
warm inflation, with just the Hubble damping term, the effective inflaton mass to^ = (y"Y^^ could be at most 
7710 ~ 10^^^" GeV < H. The key point to appreciate here is that there are flve orders of magnitude difference here 
between the vacuum energy scale and the scale of the inflaton mass. Thus there is a huge difference in scales between 
the energy scale V^^^ driving inflation and the energy scale 777^ governing inflaton dynamics. This means, for example, 
in order to excite the inflaton fluctuations above their ground state only requires a minuscule fraction of the vacuum 
energy to be dissipated as radiation, something at a level as low as 0.001%. This gives an indication that dissipative 
effects during inflation have the possibility to pla y a noticeable role (several models of warm inflation models exploiting 
these properties exist [13, [H, H [H, [H, [H M, MMMMMMMMMMMMMMMMMM) ■ 
Of course, this energy scale assessment is only suggestive. This remains a question that only a proper dynamical 
calculation can answer. In particular, the universe is expanding rapidly during inflation, at a rate characterised by 
the Hubble parameter H, and one must determine whether the fundamental dynamics responsible for dissipation can 
occur at a rate faster than the Hubble expansion. 

Another difference between warm inflation and cold inflation is how the slow-roll conditions feed back into conditions 
on the parameters of the underlying particle model. This in particular becomes evident in the strong dissipative regime 
when T ^ 3-ff. To appreciate this point, note that in cold inflation the slow-roll conditions require inflaton mass to 
satisfy < H, and this can be a problem in realistic quantum field theory models of inflation. The reason most 
realistic models of inflation rely on supersymmetry is to help cancel quantum corrections and thus maintain the desired 
flatness of the inflaton potential. However, supersymmetry can be local as well as global. Any local, or supergravity 
theory has associated with it a Kahler potential which alters the scalar field potential [l^l . Generic Kahler potentials 
lead to inflaton masses bigger than H, thus contradicting the slow roll conditions. One must then consider only very 
special Kahler potentials that can be consistent with inflation. This restricts model building prospects, but even more 
seriously, it is very likely that these very special forms may not be produced when the Kahler potential is derived, 
as opposed to put in by hand. This is often called the "7;-problem" [H, [s^, HH [E^l- There are solutions that 
attempt to stablize the flat directions [50, 53] thus allow cold inflation. Moreover a very different proposal that can 
overcome the "77-problem" in cold inflation is D-term inflation [s^, • However all these solutions require greater 
model building details. In contrast, in warm inflation slow-roll motion only requires from Eq. Ij2.16p r] < 1 + Q, which 
means that when T > 3if the inflaton mass H <m^< (i?T)i/2, which can be much bigger than in the cold inflation 
case. This relaxation on the constraint in the inflaton mass permits much greater freedom in building realistic inflaton 
models, since the "7/-problem" is eliminated. 
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Another model building consequence differing warm inflation to cold inflation relates to the range of the scalar 
field (j) in which the inflation occurs. For cold inflation, for the simplest kinds of potentials, which also are the most 
commonly used, such a,s V = A0'*/4 and V = 171^4)^/2, calculations show that the inflaton range has to be above the 
Planck scale (j) > mp. This arises because inflation ends when cj) mp, and in order to have the desired 60 or so 
e- folds of inflation, the inflaton has to start with a value larger than mp. Although the potential can still be below 
the Planck energy density, there are likely to be difficulties from quantum gravity or supergravity effects, which are 
discussed in further detail in Sec. IVIII The upshot is that restrictions have to be placed on model building simply 
to avoid this problem. On the other hand in warm inflation, when T > 3H, the added dissipation means the period 
of slow roll necessary to obtain the desired 60 or so e-folds can be achieved with the inflaton traversing over a much 
smaller range. For example, with monomial potentials, the inflaton amplitude is below the Planck scale 4> < mp 
in warm inflation. Potentials which do not allow cold inflation can sometimes be used for warm inflationary model 
building, as will be further discussed in Sec. IVIII 



C. Fluctuations 



The sources of density fluctuations in warm inflationary models are the thermal fluctuations in the radiation 
fields. This is a substantial departure from cold inflation, where the density fluctuations arise from quantum vacuum 
fluctuations. In this subsection the inflaton fluctuations during warm inflation are related to a Langevin equation for 
the inflaton field. An intuitive argument is presented and results from the systematic derivation given. 

Inflaton fluctuations are described by the amplitude S(f>(\i, t) for comoving wave number k and cosmic time t. These 
satisfy a Langevin equation similar to Eq. (I2.13p . but now in an expanding universe [8, 11], 

(5<^(k, t) + (3iJ + T)(50(k, t) + + m^)5(t){\^, t) = C(k, t). (2.20) 

In the above equation, the noise correlator is taken to satisfy the fluctuation-dissipation relation, in which case we 
have the result 



(e(k, i)e(k', t')) = 2(3i7 + T)ra-3(27r)353(k - k')(5(t - t'). (2.21) 



The noise drives scalar flcld fluctuations with amplitude (50(k, i), defined by 

((5</.(k, t)5(t){\i! ,t)) = k~^S^'^{2TrfS^{k - k'). (2.22) 

The regime of interest here is where the zero-mode is overdamped, i.e., meaning 3H + T > m. As time progresses, 
the oscillation frequency ujk = (k^a~^ -I- m^)^/^ decreases until eventually the mode gets frozen in, similar to what 
happens in cold inflation. However, since the dissipative term in warm inflation can be much larger than that in cold 
inflation due to the T term, this freeze-out momentum scale can be much larger than that in cold inflation, which 
is ~ H. At the freeze-out time tp, when the physical wavenumber hp — k/a(tp), the mode amplitude S(j) can be 
estimated using a purely thermal spectrum, 

M\kp)^f ^^(e^--l)-^-M:. (2.23) 
Jk<kF (27r)^ uJk 27r^ 

Note that, since the solutions to the source-free equation for S(j) are heavily damped, the system looses any memory 
of the initial conditions. 

To estimate kp, one must determine when the damping rate of Eq. (|2.20p falls below the expansion rate H , which 
occurs at kp « {3H + T)H. Thus, in the strong dissipative regime Q » 1, this imphes kp ~ ^/HT. Substituting for 
kp in Eq. (j2.23|l . one finds the expression for the inflaton amplitude at freeze-out 



(2.24) 

This expression was first derived by Berera in [ll|. In the weak dissipative regime Q <^1, the freeze-out wavenumber 
kp ^ H, the latter being consistent with what occurs in cold inflation. The inflaton amplitude at freeze-out becomes 
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This expression was first found by Moss in [6| and then independently rediscovered by Berera and Fang in In 



both cases it was incorrectly asserted to be the expression for the entire dissipative regime, and in the appropriate 
regime of its validity, the weak dissipative regime, was clarified. 

A much more rigorous description of the fluctuation amplitude can be found in Hall, Moss and Berera [38| and Moss 
and Xiong [56J. These papers solved the Langcvin equation (|2.20p explicitly using Green's function methods. They 
also solved the full set of equations for linear fluctuations, including metric and entropy perturbations, in addition to 
the inflaton perturbations. Solving the Langevin equation gives 

S^^ ^ k-^^[{3H + r)H]^/^T, (2.26) 

at the freezout scale, in agreement with the heuristic description above. However, a new effect can be seen when 
the friction coefScient depends on the temperature of the radiation. In this case the fluctuation amplitude has an 
oscillatory dependence on scale, caused by the entropy fluctuations which are present on sub- horizon scales. 



D. Worked example 

In this Subsection an example of the quadratic potential is presented to show how inflation models are solved in 
both the cold and warm inflation dynamics. The results will also help illustrate some key features that differentiate 
the two dynamics. 



1. Cold inflation 

We start with the slow-roll equations (|2.9[) and (|2.10p with potential 

V{ct>) = ^mlc^\ (2.27) 

To solve an inflation model, one must first determine the value c/jei at which inflation ends, and then evolve the 
field backwards to find the value of the field Ne e-folds before the end of inflation, which we call 4>n^ ■ The slow-roll 
parameters Eq. (|2.1ip for this potential become 

. = ,= ^, (2.28) 

so that the slow-roll conditions are satisfied for ^ mp\/2. Inflation ends when the slow-roll parameter e = 1, after 
which the inflaton starts oscillating and the reheating phase commences. The number of e-folds is computed as 

ftEI 1 f't'EI Y 

Ne= Hdt^~-^ -dcf>. (2.29) 

Thus, for our quadratic potential we find (f)Ei — mpV2 and 

(j)l!^=2ir?p{2N^ + l). (2.30) 

Next, we choose N,, to correspond to the largest observable scales and fix the amplitude of density fluctuations to 
coincide with the observed value. The amplitude is given by 

^ (2.31) 

5 101 

where all quantities are evaluated as the perturbation exits the horizon. Using the cold inflation expression S(j) = H 
and the slow-roll approximation, it leads to 

6h = 0.52iV2 (2.32) 
mp 

Setting 5u to the observed value « 2 x 10~^ and using A''^ = 60 leads to « 6.4 x 10~^mp. 
Finally, when computing the spectral index, n^, we find the result 
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-1^2r;-6e = - — . (2.33) 

Note that in this model the inflaton background amplitude 060 > "mp and the expansion rate Hqq > m^. Both these 
features are common to such monomial cold inflation models and, as discussed in Sec. IVIIi pose model building 
problems. 



8. Warm inflation 

Consider strong dissipation T 3> with T constant, then the first slow- roll equation becomes 

T(/) + m^(/) = 0. (2.34) 
The number of e- folds of inflation between (jj^^ and (j)Ei is now given by 



I Hdt^- I (2.35) 



Inflation ends when e = 1 + Q, where Q is defined in Eq. (|2.17p . Using the slow-roll equation (|2.35p gives (jiN^ ~ 
V&{Ne + l)mpm^/T. 

Having now determined ^tv,,, we can calculate the amplitude for density perturbations at this point and normalise 
it to the observational value. For this, the same expression Eq. (|2.3ip for the amplitude is used, except now 54) is 
given by Eq. (|2.26p . Going through the calculation, with effective particle number , we find 

5h ~ O.lSA^^/Sg-i/s ' (2.36) 

Once again setting 5h to the observational value « 2 x 10^*^ and setting ~ 60 and ^ 100, leads to a normalisation 
condition, T/mp « 1.5 x 10~^. (In fact, the number of e-folds can be as low as 40 in warm inflationary models.) 
Finally, for the spectral index, we have the general expression in warm inflation (ssj 

This is slightly closer to = 1 spectrum than the value for cold inflation case. 
From these results we find that 

^« 1 ^^^^4x10^^, (2.38) 
V6 rrip T mp 

so that choosing the inflaton mass < 10~*mp will mean that the inflaton mass is bigger than the Hubble parameter, 
thus eliminating the 77-problem, and 4>qo is below the Planck scale, thus making such models amenable to particle 
physics model building. Finally note that this analysis can be extended to other potentials [3l| and to more general 
cases where 4> and T dependence is in both T and the inflaton potential [s^. 



III. THERMAL FIELD THEORY 



We turn now to the main subject of this review, which is to describe how the dissipative effects which where put 
into the inflaton dynamics in the earlier sections have been obtained from microphysical descriptions of quantum field 
systems close to thermal equilibrium. In this section we give a basic introduction to thermal field theory with the 
assumption the reader is familiar with ordinary quantum field theory. For more extensive reviews of thermal field 
theory, please see, for example, Refs. [stI, [ssI, [sqI. [60j . 
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A. Preliminaries 



In thermal field theory the quantities of interest are ensemble averages of operator expectation values. The 'tried 
and tested' approach to thermal field theory involves expressing the observable quantities in terms of propagators and 
then applying perturbation theory. This is similar in many respects to ordinary quantum field theory, except that the 
emphasis is on the evolution of operators rather than the scattering matrix. It is best in this context to regard the 
system as being always in the 'in' state, and as we shall see below this leads to a richer propagator structure than 
usual. 

We shall be using the Schwinger-Keldysh, or the closed-time path (CTP) approach [6l[ to evaluate ensemble 
averages. We take a complete sets of states ipi and V'/ along with a density matrix p. The Schwinger-Keldysh 
generating function is defined in terms of two source terms, Ji and J2, by 



Z[Ji,J2]^Y.^HpT*exp(^-i J J20^ |V/)(V^/|rexp(^i J Ji^^i;, 

'''if 



(3.1) 



where T* denotes time ordering of the operators with the smallest time on the left. 

Ensemble averages of products of (j) can be obtained by differentiation of the generating functional with respect to 
J, for example, an average like 



{T*^{xi) . . . i{Xr)T^{Xr+l) . . . ^{Xn)) , (3.2) 

with both 'time ordering' T and 'reverse time ordering' T* is obtained from r derivatives with respect to J2 and n — r 
derivatives with respect to Ji . It proves convenient to remove the minus sign in front of J2 by defining = Ji and 
P = - J2. 

Four different connected two-point functions can be obtained from the second derivatives of the generating function. 



Ga,{x,x') - -Zjj^^^^^^^ . (3.3) 



These can be placed neatly into a 2 x 2 matrix 



r (r r') = f (Tm^i^Dc {kx')i{x))c 

"'^ ' ' { {mkx'))c {T*kx)kx')). 

We recognise that Gn is the thermal analogue of the Feynman propagator Gp- The remaining combinations are the 
thermal Dyson function G22, the thermal Wightman function G21 and its transpose G12. Note that the two-point 
functions only depend on the initial density matrix, and because of this fact the formalism is sometimes called the 
'in-in' formalism. 

The propagator can be split into real and imaginary parts by introducing the real anticommutator function F and 
real spectral function p, defined by 




p{x,x') = imx),4'{x')]),, (3.5) 
F{x,x') = i({0(a:),0(x')})c. (3.6) 



The propagator matrix separates into real and imaginary parts according to 



, ( F{x,x')~ '^(j{x,x') F{x,x') + y{x,x')\ 

^''^^^'^ > ~ \ F{x,x') ~ ^p{x,x') F{x,x') + ^(t{x,x') ) ' 

where cr(a;,a;') = p(x,x')sgn(i — t') is the real and time-symmetric Wheeler-Feynman propagator. 
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B. Thermal equilibrium 

Systems in thermal equilibrium are invariant under time translation, space translation and additionally the propa- 
gators satisfy periodicity relations in imaginary time. The propagators depend on x and x' only in the combination 
X — x' , and it usually proves convenient to use the Fourier transform over space and time, replacing x — x' with (p, uj). 

The imaginary-time periodicity relations imply that the anticommutator function F and the spectral function p, in 
space-time momentum representation, are related by j60| 

F(p,o.) = -^[l + 2nH]p(p,o.), (3.8) 
where n{uj) is the thermal distribution function for inverse temperature /3, 

= • (3-9) 

This remarkable relation between F and p, means that the full thermal propagator depends only on the spectral 
functions. 

The spectral functions can be obtained perturbatively by solving the Schwinger-Dyson equation, 

- k2 - m^)Gab - ^a'Gcb = iCab , (3.10) 

where T,ab is the self-energy matrix. The tensor Cab is a diagonal matrix with entries ±1. It is used as a metric to 
raise indices a and &, and keeps track of the minus signs introduced by the reverse time ordering. 

As with the propagator, the self-energy matrix can also be expressed in terms of two functions and Sp, 

iSp = i (E21 - S12) , (3-11) 

l^F = i(S21+Si2) . (3.12) 



In thermal equilibrium, the Fourier transforms of Ep and T,p are related by a local relation just like Eq. (|3.8p . 

S;^(p, u;) = -'-[l + 2niu;)] Ep(p, u) . (3.13) 

The physical interpretation of Sp is that it is related to the decay processes and can be associated with a relaxation 
time defined by 

r(p,c.) = — (3.14) 
«Sp(p,tj) 

Equation (j3.13p will turn out later to be related to a fluctuation-dissipation theorem. 

The real time formalism extends to Dirac spinors in a routine manner. A description of the thermal propagator for 
Dirac fields can be found in [605 ■ The propagator in similar conventions to the ones used here can be found in [l7j . 

C. Real time formalism for interacting field theories 

We now explain how to relate the general Green's functions to the free field case using perturbation theory, i.e., 
Feynman diagram expansion, in the context of the Schwinger-Keldysh formalism. The causality properties of the 
Feynman diagrams are also discussed. 

In path integral form, the generating function p.ip becomes 

Z[Ji,J2] = J dfi[(l)i]dfi[cl)2]p[MU),MU)] exp|z5[(/.i]+i J Ml ~ iS[<j>2] - t J J2hy (3.15) 
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where p is the density matrix at the initial time ti and the paths cross asymptotically as t — > oo. This path integral 
is equivalent to using a single field on a closed time path (CTP), taking the time integration along a contour in the 
complex time plane going from t — ti to +00 (forward branch) and then back io t — ti (backwards branch). We 
still refer to the formalism as the CTP formalism even though we find it convenient to keep the fields and branches 
distinct. 

The CTP Feynman diagram expansion for ensemble averages mirrors the ordinary Feynman diagram expansion 
for n-point functions very closely. An important difference is that vertices carry an extra integer label taking the 
value 1 or 2 which determines the component of the propagator matrix between them. The vertices represent the 
difference of interaction Lagrangians, Ci{4>i) — J~-ii4'2), and therefore a minus sign is included for each vertex labeled 
2. Troublesome minus signs can be removed from the source terms by using a metric — diag(l, —1). 

One useful feature of the Feynman diagram expansion is the maximum time rule. If the diagrams are drawn in 
configuration space, then a diagram gives a zero contribution if the time on any internal vertex is larger than the 
largest external time. This rule follows from the symmetries of the Green's function. It implies that the effects of 
source terms always satisfy the rules of causality, event though the Green's function is acausal. 

The CTP effective action also follows by direct analogy with the usual effective action, 

e^n4>i,M^ f dA*[0'i]dAi[(/)'2]p[0'i(i.),(/.'2(t.)] exp{z^[0l + 0i]-i5[0^ + 02]}. (3.16) 
JiPi 

where only the 1-particle irreducible (IPI) diagrams contribute. It satisfies the effective field equation 

& =0. (3.17) 

The new features introduced in the CTP approach are the doubling of fields and the condition 0i = (f>2- 

The path integral may be constructed from the original vertices and Green's functions, or by shifting the Lagrangian 
about a background field and then using shifted vertices and corrected Green's functions. When the propagators and 
the vertices of the Feynman diagram expansion depend on the background fields, they must be kept distinct. However, 
after taking the variation of the effective action to obtain the effective field equations, the values of (jji and 4>2 are set 
to the same value. Note that the maximum time rule quoted above only applies when (jji = 4>2- 

An important case of the background field approach is when the system is close to thermal equilibrium and the 
background fields vary slowly compared to the relaxation time of the system. In this case, it can be appropriate to take 
the propagator of the background field expansion to be in thermal equilibrium, and relate all of the non-equilibrium 
effects to the background field dependent interaction terms. 



IV. THE EFFECTIVE EQUATIONS OF MOTION 



Now we introduce effective equations of motion, where there is a background field, which makes the role of a system, 
and quantum fields, considered as the environment to which the system is coupled to. By integrating out the quantum 
fields it is possible to arrive at a background field equation, which is of a Langevin-like type typical of a system in 
interaction with an environment. We will keep our formalism as simple and general as possible, allowing then to 
extend it to specific model examples later on. 



A. Historical background and motivation 



The basic motivation for the study of the nonequilibrium dynamics of a background field in the context of a 
separation between a system and an environment to which it is coupled to comes from many different physical 
systems of interest. In the case of inflation, it is the inflaton field whose slow-roll dynamics generates the necessary 
conditions for a successful inflationary scenario. Almost all models of inflation involve evolution with loss of energy of 
the inflaton field to other fields (or degrees of freedom) to which it is coupled to. This process of energy transfer from 
the inflaton (regarded as the system) to the other fields is a fundamental requirement of equipartition, where some 
portion of the system's energy has to fiow irreversibly to the environment. This process is not exclusive to infiation, 
and occurs in any phase transition where some order parameter characterising the global thermodynamic properties 
of the system relaxes to an equilibrium point. 

The study of the nonequilibrium dynamics of fields have been approached using different techniques of quantum field 
theory and quantum statistical mechanics, including variational techniques [6^ . the use of resummation techni ques , 
for example using the two-particle irreducible (2PI) procedure ^GSj], and the use of kinetic equation methods [64| . 
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These different approaches have in common the possibihty of keeping all the dynamics unitary, but at the expense of 
keeping track of every field, and not just the one whose dynamics we are most interested in. Even though this study 
might be done in a relatively complete fashion for some simple models, it gets quickly cumbersome as the number of 
fields (and field modes) increases. 

In the Langcvin-likc approach we focus on the dynamics of the relevant field describing the system and not on the 
the remaining field modes. These instead act on the system through dissipative and stochastic noise terms, whose 
effects then become manifest. This is a much more economical way of studying the nonequilibrium dynamics, since 
we concentrate only on a given (more relevant) field mode. This approach has a long history which can be traced 
back to Einstein's famous explanation of Brownian motion [6^, and includes classics such as the work of Caldeira 
and Leggett [6^. In their work an harmonic oscillator, regarded as the system, is coupled linearly to a set of other 
oscillators, regarded as the thermal bath or environment. By integrating out the thermal bath, the resulting dynamics 
of the system becomes explicitly of Langevin form, where dissipation and stochastic noise emerges. Extension of these 
studies to the context of nonlinear couplings between system and thermal bath have been implemented by Hu, Paz 
and Zhang [67l |. 

Though in quantum mechanics the system, which is out of equilibrium, and the thermal bath, which drives the 
system towards equilibrium, may be well separated, in the context of nonlinear field theories this distinction may be 
considered somewhat blurred. Even so, for self-interacting field theories, there are situations where short wavelength 
modes can serve as the thermal bath driving the longer wavelength modes, which have slower dynamics, into equilib- 
rium. In this sense, the field can be its own thermal bath. Of course, other fields coupled to a background scalar field 
(the system) may also serve as the thermal bath. One of the first implementations of this interpretation in the context 
of quantum field theory, and motivated by the reheating problem in inflation, was the work by Hosoya and Sakagami 
[6^, who obtained an approximate dissipation term in the equation of motion for a scalar field. They did this by 
examining small deviations from equilibrium in the Boltzmann equation for the number density operator and then 
supplemented this by a computation of transport coefficients using Zubarev's method for nonequilibrium statistical 
operators [HH] . The derivation of dissipation terms in the context of the A(/>^ model was also performedusing operator 
methods by Morikawa and Sasaki in 69]. Later, in the context of the OTP formalism, Morikawa 70*| obtained an 
effective Langevin-like equation for a scalar field interacting with a fcrmionic bath, including explicit fluctuation and 
dissipation terms. 

Other work that analysed the emergence of dissipation and fluctuation in the effective field dynamics was the work 
done by Hu and collaborators (Tlj , who analysed a scalar field quantum bath quadratically coupled to a background 
scalar field system, while Lee and Boyanovsky considered the case of a scalar field thermal bath linearly coupled to 
a background scalar field system [t^, with more realistic couplings considered later on [t^. In the work of Gleiser 



and Ramos |74| , a systematic study in the context of the loop expansion at high temperature was performed for both 
the model and also for quadratic coupling to another scalar field. One important aspect of the dynamics that 
was demonstrated by these first references was that the noise terms emerging in the effective dynamics was in general 
colored, i.e. non Markovian, and multiplicative, i.e. field dependent (unless the coupling between system and bath 
was linear). Thus, these studies have shown that the effective equation that describes the approach to equilibrium of 
the slower moving modes can be quite different from the typical phenomenological Langevin equation with its white 
and additive noise terms. 

All these early studies were performed in the context of Minkowski spacetime. One of the first to consider dissipative 
dynamics in the context of a curved spacetime was Ringwald [TSj] . Later on, the problem of dissipation and damping 
in a de Sitter spacetime was considered by the authors of Ref. [7y| , while more recently the nonequilibrium dynamics 
of the infiaton field in the Friedmann- Robertson- Walker spacetime was considered by Berera and Ramos in [73] , where 
an extensive analysis of the dissipation kernels, entropy and particle production were performed. In this section and 
in the next one we will be most concerned with the dynamics in Minkowski spacetime, while in Sec. I VII we will discuss 
the changes necessary to implement in order to describe our results in the context of the Friedmann- Robertson- Walker 
curved spacetime. 



B. Introducing the Keldysh Representation 

We start by noting, from the definitions given in Sec. IIIII for the four two-point functions defined in the CTP 
formalism, Eqs. (|3.7p . that they are not all independent. This indicates that we can define a linear transformation of 
the fields to make some components of the propagator matrix vanish. In terms of the field defined in the forward 
branch of the CTP contour <j)i and in the backwards branch (j)2, this linear transformation (also called Keldysh 
rotation) leads to two new fields (j)c and (pA defined by 
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+ 02) 



[>2 ■ 



Now, it is a matter of simple algebra to show that the propagator matrix transforms to (a', b' — c, A) 



(4.1) 
(4.2) 



^ _ / F{x,x') Gr{x,x') 
- \Ga{x,x') 

where 



(4.3) 



Gr{x, x') = -ip{x, x')0{t - t') 
Ga{x,x') = ip{x,x')e{t' - t) , 





'iln{x,x') 



Y.Aix.x') 
-iT,F{x, x') 



where 



= j:p{x,xy{t-t'), 
Ea = -j:p{x,x')e{t' -t). 



(4.4) 
(4.5) 



Gr and Ga can be identified with the retarded and advanced two-point functions, while the function F{x,x'), defined 
in Eq. (|3.6p . is sometimes also called the Keldysh two-point function. 

The self-energy matrix in the Keldysh representation also has a simple form. 



(4.6) 



(4.7) 
(4.8) 



Let us see some of the advantages of working with the Keldysh representation for the fields. As an example consider 
the classical action for a theory in the CTP formalism, 



d^X 



d^x 



-01 (-92 _ ^2) 



A 

4!' 
A 



'2 - ^<P2 



(4.9) 



which in terms of (/>c, (pA becomes 



OA 



d'^X 



-92 



A 
4! 



(4.10) 



The first thing to note from the action (/ia] is that it vanishes for a field configuration that it is the same on 
the forward and backwards branches of the closed-time path, i.e. when <f>\ = 0. Although this seems obvious for the 
classical action, it is important to realise that the same structure remains true for the effective action, with no terms 
independent of 0a appearing at any perturbative order. This restricts the form of vertex and self-energy corrections. 
For example, the c — c term in the self energy has to vanish because of the form of the self-energy matrix, Eq. ()4.6|) . 
Another consequence of using the representation of the action in terms of the (pc and (/)a fields is that the classical 
field equation can be obtained from 



0A=O 



= 0, 



(4.11) 



which for Eq. (|4.10p can be seen to immediately reproduce the usual classical equation of motion for a A(/)*/4! theory. 



(92+m2)0, + A03^O 



(4.12) 



The down-side of working with the Keldysh representation is that the simplicity of the propagators is offset by the 
complexity of the vertices in the Feynman diagram expansion. The best approach often is to use the original matrix 
propagator of the Feynman diagram expansion and then to use the Keldysh representation for the end result. 
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C. The Effective Action and Equation of Motion: A System-Environment (Langevin) Interpretation 

Due to the special properties of the Keldysh representation, the quantum corrections to the classical action, i.e. 
the effective action Eq. (|3.16p . can be represented generically in the form: 

r[(bc,M = - J d^xJ^{x)<PA{x) + ^ J d^xd^x' <j)A(.x)ii:F{x,x')q^Aix') + 0{q^l) , (4.13) 

where we have kept terms up to second order in the field 0a • As we are going to see next, these terms have an 
important interpretation in the definition of the effective dynamics for a background field. 

The first term in Eq. (|4.13p can be recognized as the term which leads to the effective field equation T ~ Q, in the 
absence of the second term in that equation, where 



nx) = 



ST[(l)c, (pi 



OA 



(4.14) 



f>&=a 



Perturbatively, this is given by the classical field equations plus corrections from the 1-particle irreducible Feyn- 
man diagrams. For the moment, we shall drop the vertex corrections but keep the full the self-energy. (A better 
approximation scheme is adopted in the next section.) The effective field equation becomes 



3! 



(blix) 



Mx)+ / d^x'^Bix,x')(f,c{x')=0 



(4.15) 



where the self-energy term Sfl(a;, x') = T,p{x, x')6{t — t'). An important property of the effective equations of motion 
in the Keldysh formalism is that the causality is always explicit, a fact reflected here in the use of the retarded 
combination of self-energy. 

The second term in Eq. (|4.13p is a purely imaginary term in the effective action that depends only on the self- 
energy. This term contains information which is needed to describe fluctuations about the solutions to the effective 
field equations. A useful trick has been developed which replaces the quantum fluctuations by statistical fluctuations 
in an ensemble of random fields. This is done by performing a Hubbard-Stratonovich transformation in the functional 
partition function, introducing a random field ^{x) to decouple the quadratic term in (f>A in Eq. (|4.13p . 

Consider a functional integral with the classical action replaced by the effective action (|4.13p . The tree diagram 
contributions to this functional integral generates the full n-point functions, just as it does in the non-CPT formalism. 
The quadratic term in (/)a in the functional integrand can be written as 



exp| — - J d'^xd'^x' (f>A{x)T,p[(t)c]{x, x')(f>A{x') 

= |detSj.|i/2y exp|-i J d'^xd^x' ^{x)Y.p\x,x')^{x') + t J d'^x^{x)(f>A{x)^ , (4.16) 

The functional integrand now has a real quadratic term in the field ^{x), and a linear term in (j}A- Taking the tree 
diagram contributions to the functional integral with this new effective action gives a stochastic equation of motion, 



d'+m' + ^_Plix) 



Mx) + / d^x' 1:r{x, x') 4>c{x') = i{x) . 



(4.17) 



Equation (|4T7|) can be seen as a Langevin-like equation of motion. From Eq. (|4.16[) . ^(x) can be interpreted as 
a Gaussian stochastic noise with the general properties of having zero mean, {^{x)) ~ 0, and two-point statistical 
correlation function 



{axmx'))^^F{x,x'). (4.18) 

Statistical averages are defined as functional integrals over the ^{x) field. An important property of Eq. (|4.17p . related 
to the Langevin-like form, is the existence of a dissipative-like term. To demonstrate this, note that we can define a 
dissipation kernel I?(a;, x') as [tI, [T^I 
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dt 

and Eq. (|4.17|) then becomes, 



^pix,x') = -^Vix,x') , (4.19) 



(4.20) 



Under a space-time Fourier transform, using the definition (|4.19p . we find that the noise kernel given by Eq. (|4.18p 
and the dissipation kernel I?(a;, x') in Eq. (|4.20p are related by the relation (|3.13p . 



Ei?(p, uj) — 2u! 



n{uj) + - 



(4.21) 



Note that, in the Rayleigh- Jeans regime < T, 2w [n(w) + 1/2] ^ 2T, so that Eq. (|4.2ip and Eq. (|4.18p reproduce 
the classical relation between the fiuctuation two-point function and the dissipation. 



(^(P, t')) = 2T 1 ^Vip, c.)e'"(*-*') , (4.22) 

which forms the basis of Eq. (|2.2ip used in Sec. II. In Ref. was demonstrated that, at high temperatures 

(typically T ^ m^, w^), the noise and dissipation kernels tend to approach local forms and then Eq. (|4.20p becomes 
Markovian, with a white multiplicative noise term. Improved analysis for the localization of the noise, and consequently 
the dissipation term, was recently done in Ref. (79j . 



V. MARKOVIAN APPROXIMATION FOR THE NOISE AND DISSIPATION KERNELS 

The possibility of approximating nonlinear and nonlocal equations of motion like Eq. (|4.20p in a local form offers 
many advantages. Typically, solving nonlinear stochastic equations of the form of Eq. (|4.20p is very hard both 
analytically and numerically. Any method attempting to solve these kind of equations requires keeping the memory 
of the past history of the scalar field configuration at each stage of the evolution. These equations also typically 
involve highly oscillating nonlocal kernels that can lead to errors which quickly build up and that are too hard to 
control, thus preventing any simple numerical solution. There is an immense saving of effort as well as a much better 
understanding of the physics from a local equation as opposed to a nonlocal one, since the former can generally be 
analysed with a more transparent numerical treatment than the latter. For these reasons, attempts have been made 
to express the equations in a local or Markovian approximate form. In particular, it has been suggested in [isj that 
at high temperatures and with a large set of heat-bath fields, the existence of many decay channels could lead to an 
approximate local Langevin equation of motion for (f>. The results in [Tsj motivated one of the first microscopically 
motivated models for warm inflation [l^. The large set of heat-bath fields proposed in [Tsj would constitute of a 
tower of massive modes, in a string motivated model, through which the inflaton field could interact. In the model 
proposed in enough radiation would be produced, leading to an overdamped motion for the inflaton and making 
possible to sustain inflation long enough. 

Later attempts to treat the nonlocal kernels, and not relying on a local approximation, but still at high temper- 
ature with a large thermal bath, where proposed in Refs. [lg|. There the analysis were based on the possibility of 
constructing models where the kernels exhibit a strong exponential damping in time, making possible a numerical 
approach to analyse the nonlocal equation of motion for the scalar field background. The results in Refs. [l^ have 
also shown that a local approximation for the kernels, and then for the equation of motion, is in very good agreement 
with the full numerical solution of the nonlocal equation. However, as shown in Ref. (Tt} . there are cases where such 
a strong damping behavior for the kernels are not possible and, thus, we must resort to other alternative analysis to 
determine how good a local approximation is for the nonlocal equation. 

Local approximations for equations of motion of the form of Eq. (j4.20p have been criticized by a number of authors. 
Lawrie [80|, for example, has studied the non-equilibrium dynamics using various approximations to the propagators. 
He has argued, from a formal point of view, that the local approximation would violate some specific sum rules in 
the kinetic equation approach. In the cases where the local approximation was tested against the numerical solutions 
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coming from the kinetic equations derived in [80( 1 , it was found that the local approximation tended to over-estimate 
the real dissipation. A similar conclusion was reached by Aarts and Tranberg [sj], who used a numerical code to 
evolve the propagators in a large W approximation for various models which resemble those used in warm inflation. 

The main drawback with these numerical approaches so far is that they have, of necessity, to be based on models 
in which warm inflation is not expected to occur even in the closc-to equilibrium approximation. This is because they 
have no mechanism to suppress thermal corrections to the inflaton potential and, as we arg ued in Sect IIIB2[ warm 
inflation cannot take place. The warm inflation models discussed in Sec. I VIII all use some form of supersymmetry, and 
at the present time no fully non-equilibrium calculation has been possible due to the complexity of the field content. 

Very recently, in Ref. [79| , specific conditions were derived for the validity of adopting a local approximation for the 
dissipation and fluctuation kernels in a specific model favored by warm inflation and that can be physically realised 
in the context of supcrsymmctric models [TtI . [s^ . 83] . In this model, first identified in [T6| , the background scalar 
field (j)c is coupled to heavy intermediate quantum fields which in turn are coupled to the light quantum fields. The 
dynamics of dissipation and radiation production in this model is realised by a two-stage mechanism: the background 
scalar field indirectly induces particle production in the light fields through the intermediate heavy fields which in 
a sense help to catalyze the effect. A throughout study of model realisations of such mechanism for dissipation, in 
systems near thermal equilibrium has been given in Ref. [l7j , while application of these effects to inflation have been 
shown to have significant importance [4^ . 

One important result drawn from the analysis of Ref. |79j was the demonstration of the nonexistence of a local 
approximation for the dissipation and fluctuation kernels at zero temperature. This result was shown to be a direct 
consequence of the application of the Markovian approximation to the generalised fluctuation-dissipation relation like 
Eq. (|4.2ip . This in particular implies that there should be no local first order time derivatives in the equation of 
motion like Eq. (|4.20p at T = 0. This result was also explicitly shown to be the case in Ref. [l3| (though it was also 
shown that higher order local but non analytic derivative terms were possible). 



A. Dissipative Effects: Local Approximation 

As discussed above, the dissipative term takes a local form when the background field is slowly varying and the 
system remains close to thermal equilibrium. The dissipation in this case is related to the transport coefficient which 
we have been calling the friction coefficient T. We now give a general formula for the friction coefRcient. 

Consider the effective field equation for the field in the CTP formalism, which was generated by 



Ha{x) 



(5.1) 

(«A=0 



In the following, let us explicitly consider that (t>c is spatially homogeneous and that it varies slowly about its value 
(/)(i) at a fixed time t. Set 5(j)c = 4>c — 4'{i) expand T by 



oo 

.F(x) = ^.F„(x), (5.2) 

n=0 

where 



1 



d^Xi . . . d^Xn 



'>Aix)S(t)c{xi) . . . 6(t)ciXn) 



(5.3) 



The first term J^q represents the part of the field equations which contains no derivative terms, and can be expressed 
as the derivative of an effective potential V{4>), 



54)a{x) 



dV 



(5.4) 



The next term Ti depends on the equilibrium self-energy of the inflaton with constant values of the background 
field, 



5(t)^^{x)5(j)c{x') 



= {d^ + m^)S{x - x') + Y.ii{x, x') , 

0A=O 



(5.5) 
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where I]_r = T,p6{t — t'), as before. The total field equation up to first order in S4> becomes 

4> + J d^xi Enix - xi)S(f>ih) + 1^ = 0. (5.6) 

We re-iterate that, when the self-energy is calculated, we can take (fi to be constant. Since we expand about thermal 
equilibrium, we have used the fact that S/j(x,Xi) = 'Sj^^x — Xi). 

The non-local dissipative term can be localized when there is a separation of timescales in the system. Suppose, 
for example, that the self-energy introduces a response timescale r. If is slowly varying on the response timescale 
T, then we can use a simple Taylor expansion and write 

= (/.(t) + (ti-i)0(t) + ... (5.7) 

The (f> equation of motion including the linear dissipative terms is then 

■ dV 

(/)-f T(^-f— =0, (5.8) 



with dissipation coefficient 



(5.9) 



where Sp is related to Tj^. in Eq. (|4.7p . This equation forms the starting point for calculating the friction coefficient 
in particular particle models. The general expression relates the friction coefficient to the imaginary part of the self 
energy Ep(p, w) at zero momentum. 



B. Coupled field systems 

We shall now obtain the friction term in a basic example with an inflaton field, another scalar field X; and thermal 
radiation fields a. The field x in this example acts as the only channel for the transfer of energy from the inflaton 
into heat radiation. This situation offers the best prospect so far for realising warm inflation in realistic models, as 
it gives some degree of separation between the processes which govern the thcrmalisation of the heat bath and the 
coupling constants which are relevant to the dissipative dynamics of the scalar field. Model building will be discussed 
further in Sec. IVlIl 

A suitable Lagrangian for the infiaton and x interactions is 

Ci=gmm + SrM'' + 2g^m''\x\\ (5.10) 

where we use </> to denote the background value of the inflaton field and 5(j) to denote the fluctuating components. 
For the x a-nd cr interactions. 



C'^^—hm{o\*+a*\). (5.11) 

Complex fields are used because they embed more easily into supersymmetric theories, which, as discussed previously 
in Sec. |lTl are more suitable to describe realistic warm infiation models by keeping quantum (and thermal) corrections 
to the potential, that would be otherwise harmful, small enough. 

The contribution to the self-energy of the inflaton field at order is given by the first diagram in Fig. [3] with two 
X propagators. When the self-energy is expressed in terms of the spectral function of the x field, one obtains a 
formula for the dissipation coefficient. 



For small h and fixed T, the energy integral is dominated by the point ojk = (fc^ + m^)^/^, which lies close to two 
poles in the spectral function. These two poles are at w = cj^ ± *t~^, where is the relaxation time for the x boson 
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FIG. 3: Contributions to the (f> self-energy of order (left) and (right). 

defined by Eq. (I3.14p . The integrand can be expanded about u — ojk to obtain a formula first obtained by Hosoya 
and Sakagami [68,], 

T«ffW/3|^^n(n + l). (5.13) 

The dependence on a relaxation timescale is typical of the structure which one might expect from elementary transport 
theory. Note that reducing the coupling constant h increases the relaxation time (since ~ 0{h~^)) and therefore 
increases the friction coefficient. This is also typical of elementary transport theory and is seen, for example, in the 
Drude theory of conductivity. The obvious 'reductio ad absurdam' argument of reducing the coupling to zero does not 
apply because our assumption that the system remains near thermal equilibrium sets an upper limit to the relaxation 
time. 

The contribution to due to the single a loop with coupling h gives ujt~^ = h^ni^ /{32tt), and [isl. [68l. [t^ 

T^—^Tln—, T:^m^. (5.14) 

Note that including other interactions which reduce will also reduce T. 

The approximation used to derive Eq. (|5.13p fails at low temperatures, when the low energy and momentum 
behavior of the spectral function becomes the crucial consideration. This case was first addressed correctly in [l^ . In 
cases where the x-field self-energy is non- vanishing in a neighborhood of A: = w = 0, we use the relation « (T,p)^/m^ 
to deduce that 



T^Cg'h*(—) ^, T«m^, (5.15) 

for a constant C, which can be determined accurately from numerical integration. For the interaction Lagrangian 
given above, C ~ 0.006. In figure IH we plot the overall behaviour for the dissipation coefficient for different regimes 
of temperature. 



A fully supersymmetric theory has many other interaction terms which contribute to the dissipation term in the 
inflaton equation of motion. At high temperatures we must consider an interaction term. 



Ci = 2h'\x\'\<jf, (5.16) 

in addition to those in Eqs. (|5.10p and (|5.1ip . This term dominates the X'fi^ld self-energy at large temperatures 
leading to ujt~^ ~ /i*T^/(1287r^) and, from an analogous expression to the one given by Eq. (|5.13p . it gives for the 
dissipation coefficient the result [H, [13, [ill 



T l-m^T-Mn , /iT > m^. (5.17) 
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FIG. 4: The difFerent approximations for the friction coefficient are shown in the intermediate and low temperature region. 
The full expression, plotted on the left, corresponds to Eq. (|5.12[) and the low temperature approximation to Eq. (|5.15p . These 



plots include both 



2cr and 4> decay channels. Coupling constants are h /8tv = 0.025 and 



Another interaction which may be present in the supersymmetric theory is 



Ci^gh{dc^X + 5rx*)W\^- (5.18) 

This allows a direct interaction between the inflaton and the heat bath, but it does so without affecting the infiaton 
potential in any serious way. The interaction results in a vertex correction to the inflaton self-energy which modifies 
the friction coefficient T, increasing the value of the constant C, in Eq. (|5.15p . to C « 0.023. 

Fermionic decays can also be included, replacing the scalar field x by a fermion or replacing the scalar field cr by a 
fermion. These are also considered in [I3, where the different dissipation coefficients for each case can be found. The 
fermionic heat bath fields reduce the high temperature friction coefficient, but at low temperatures they contribute 
terms (x to the friction coefficient, which are negligible compared to the bosonic contributions. 



C. Physical picture of particle production 

So far we have concentrated on the dissipative effects produced by the interaction between the inflaton and the 
radiation fields. The underlying process here is some form of particle production, and it should be possible to see the 
same physics by a consideration of the particle production rates. This is in fact the case, and this approach offers 
a way of understanding the thermalisation processes in the radiation, or indeed, of discussing what happens when 
thermalisation is incomplete. 

Energy conservation implies that the evolution of the total radiation energy density is given by 



= - AHpr. (5.19) 

This equation may have an equilibrium point where the redshift and the particle production vanish. In the thermal 
case, this equation has a stable equilibrium when T oc T'^, but not if T oc or any higher power of the temperature. 
It was therefore very important that the friction coefficients given in the previous subsection had the necessary 
behaviour. 
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We can weaken the thermal assumption and assume a quasiparticle approximation in which the propagators have 
a similar form to thermal propagators but where the momentum distribution function of the radiation fields n{p,t) 
is arbitrary. The energy density is then 



Pr = J J2^n{p,t)ujp. (5.20) 

The distribution function for the radiation evolves by a Boltzmann-type of equation with a source term representing 
particle production Sp from the evolving inflaton fields, a collision term Sc due to the field interactions and a redshift 
term Sr caused by the expansion of the universe, 



n, = Sp + SR + Sc ■ (5.21) 

Just as we discussed for the total energy density, equilibrium may occur when the redshift and the particle production 
balance, but now we also need the collision integral to drive the momentum dependence towards a thermal spectrum. 
Given a sufficiently large self-coupling for the radiation fields there is no reason in principle why thermalisation 
cannot occur, and numerical solutions support this conclusion [83| . An interesting possibility is that departures from 
a thermal distribution can be studied using this approach and their effects on the density fluctuations analysed. This 
problem is also well-suited to numerical analysis. 

The source term depends on the details of the particle production mechanism. For the two-stage decay mechanism 
used in Sec. IV B( the slowly evolving inflaton field cannot produce very massive x particles directly, but it can decay 
into massless radiation fields via an intermediate virtual x channel. The source term for massless radiation can be 
found analytically (841] . 



= ^F(p)0^ T«m^, (5.22) 

where F(p) is plotted in figure[5] The momentum distribution is larger at low momentum than a thermal distribution 
with similar total energy. Integrating the Boltzmann equation over momentum recovers the energy density equation 
(|5.19p with the same friction coefficient as was obtained before in Eq. (|5.15p . This provides an important consistency 
check between the dissipation and the particle production mechanisms. 



VI. EXTENDING DISSIPATIVE DYNAMICS TO CURVED FRIEDMANN-ROBERTSON- WALKER 

SPACE-TIME 



In order to complete our discussion of dissipative dynamics we now consider how the results quoted so far extend to 
curved space-time. Specifically, we consider a homogeneous and isotropic, spatially flat, Friedmann-Robertson- Walker 
metric with scale factor a. The naive expectation would be that thermal effects are more important than curved space 
quantum effects when we have radiation with temperature T which is much larger than the expansion rate H,T ^ H. 
This expectation can be supported by an approximation scheme given below, and therefore curved space quantum 
effects can be considered small during warm inflation, at least during a thermal regime. 

Scalar field propagators can be defined as in Eq. p.4p on the Friedmann-Robertson- Walker background and satisfy 
a curved space version of the propagator equation. 



[dj + 3Hdt - a-^V^ + m^+ (Rit)] Gabix, x') + j d^y Sa'(x, y)Gcb{y, x) = iCatSix, x') , (6.1) 

where R is the curvature scalar, R = 6H + 12H^, and C is dimensionless parameter describing the coupling of matter 
fields to the gravitational background. Both d'^y and S{x, x') implicitly include \/—g factors, where g is the metric 
determinant, e.g. 



5{x, x') — 



5\x - x') 



(6.2) 
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FIG. 5: This plot shows the momentum dependence of the particle production rate Sp for the production of low mass fields 
through an intermediate heavy field. The thermal distribution n is shown for comparison. 

as required by general covariance. 

The Lagrangian density for the inflaton field given in Eq. (I4.9[) is now, 



m = 



2 



4! 



(6.3) 



We can decompose the propagator and self-energy matrices as we did in flat space. Thus, for the Lagrangian density 
model given by Eq. (|6.3[) . it follows by the same arguments given in Sec. IIVI that the effective stochastic equation of 
motion for the background scalar field (pc is given by 



9? 



A 



(6.4) 



where noise field ^(x) again satisfies {£,{x)) = and it has the same two-point function as in Eq. (|4.18p . but where 
the self-energy terms are expressed in terms of propagators in curved space-time. 

These expressions hide the fact that the propagators for the fields in curved space-time are in general very com- 
plicated functions [H, [sH]. For the situation of interest here, we have a homogeneous background with scalar and 
radiation fields present. In the spatially homogeneous case, let 



Gafc(k, t, t') - j d^X Gab{x, X') e^k (x-x') (g 5) 

We shall make a pseudoparticle (or Kadanoff-Baym) approximation for the radiation (and x field) propagators, 
introducing an occupation number n(k, t). This may, for example, be thermal at some initial time to with frequency 
Lulto), where 



(6.6) 



with M the quadratic mass term for the field in the curved space-time. The spectral and anticommutator functions 
defined by Eqs. p.Sp and p.6p . respectively, in the pseudo-particle approximation are 
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F(k,t,i') = 



[/i(k,i)/2(k,i')-/2(k,t)/i(k,t')] 



n(k,t) + - 



[/l(k,i)/2(k,t') + /2(k,i)/l(k,i')] , 



(6.7) 
(6.8) 



where the functions /i,2(k, t) are defined by the solutions of the differential equation [87 



d „d d , ,9/ X 



/i.2(k,<) = 



(6.9) 



with Wronskian /i(i)/2(i) - /i(0/2(0 = -i/a'(t). 

Usually, solutions for Eq. (|6.9p are known only for some S2ecific cases, e.g. for de Sitter expansion H ^ constant, 
so a{t) = exp{Ht) [sl], and power law expansion a{t) ~ [s^. For example, in de Sitter the solutions for constant 
M are given in terms of Hankel functions, 



/i(k,t) = mk,t) = ^i7-V2e-3H*/2^(i) ^ke-"'/H) 



(6.10) 



with = 9/4 - M^/H^ - 12C0. 

Alternatively, adiabatic approximations for the mode functions can be derived by applying a WKB approximation 
to Eq. (EH) [si, [13, 



/i(k,t) = /*(k,i) 



a3/2(t)^2^I) 



exp 



dt'ioit') 



(6.11) 



Note that this leads to a well defined split between positive and negative frequency modes. 

The WKB approximation requires the adiabatic conditions ^ w^/w^, w/w, which apply at large mass or large 
momentum. Both of these are relevant to models of warm inflation, where we may have large mass (for the x held) 
and large momentum due to large temperature (compared to the expansion rate), as shown in Sec. IVIII 

These adiabatic conditions also_simplify the equations for the propagators and allow interactions to be taken into 
account. Berera and Ramos in [77] have demonstrated both numerically and also analytically that the dynamics of 
the inflaton including dissipation effects, is well approximated by the Minskowkian dynamics, thus allowing us to use 
the results in Subsec. IV Bl for the dissipation coefficients in warm inflation. 

The energy density of the radiation pr in the pseudo-particle approximation can be expressed in terms of the 
anti-commutator function (|6.8p . For massless radiation in de Sitter space, the energy density is given by 



Pr — Pdc Sitter + 



d^k 

(2^ 



1 



2i72/c2 



k n(k, t) , 



(6.12) 



where pdc Sitter is the energy density of the de Sitter vacuum. A consequence of this expression is that the energy 
density rapidly redshifts towards the de Sitter value if n(fc, t) is time independent. On the other hand, the combination 
of particle production and particle interactions can produce a thermal distribution in the physical energy spectrum 
as in Sec. IV CI i.e., n = n(k/a) and then we recover a thermal contribution to the energy density. 



VII. PARTICLE PHYSICS MODEL BUILDING FOR WARM INFLATION 

Several phenomenological warm inflation models have been constructed [13, IH, [H, [13, [H, [1^ [13, [lH, [13, [H, [H, 
[H, [H, [13,IH, [iOj [U [IE [IE [ill ■ Moreover, interesting applications of the warm inflation regime have been suggested 
for generating cosmic magnetic fields [26*1 and for baryogenesis [13| . More interesting are models of warm inflation 
constructed completely from first principles quantum field theory. There are by now many such models and these 
models have some unique and attractive features. First, in the strong dissipative regime they have no "?7"-problem. 
The "ry"-problem typically emerges since supergravity (SUGRA) corrections to the inflaton mass are of order the 
Hubble scale, yet in cold inflation the inflaton needs a mass less than the Hubble scale to realise slow- roll inflation 
[isl. [4^ [sol. [Sll. Is^ . In contrast, in strong dissipative warm inflation, the inflaton mass is much bigger than the Hubble 
scale, so such models are fairly insensitive to SUGRA corrections. The second attractive feature of warm inflation 
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models is for monomial potentials the amplitude (j) of the inflaton field is always below the Planck scale. In cold 
inflation, for monomial potentials the inflaton amplitude during inflation is larger than the Planck scale. This is 
a problem for model building. Quantum field theory models are generally regarded as low energy effective theories 
of some higher more fundamental theory, such as possibly strings, and there is some upper energy scale to which this 
low energy approximation is valid. Above this scale, the theory would be modified by additional, usually an infinite 
number of, operator corrections. The highest scale can be the Planck scale, and so for any quantum theory above 
this scale one expects an infinite number of nonrenormalisable operator corrections, ^ X^J^Li ff"'/'^(<^/'^J')": which 
have to be retained p^ . Isol [sil . [s^ . In such a regime the low energy approximation to the theory is essentially not 
useful. Thus from a model building perspective cold monomial inflation models are difficult to implement, whereas 
warm inflation models discussed in this section do not suffer this complication. 

Hilltop inflation models [90, 91} are one type of single field cold inflation model where the inflaton amplitude remains 
below the Planck scale initially during inflation, thus when the large scale features are determined. However even 
in these models typically the inflaton amplitude goes above the Planck scale by the end of inflation. Beyond these 
models, maintaining the inflaton amplitude below the Planck scale in cold inflation models requires more elaborate 
constructions, with a few examples such as [H, [sl, [9^1 • However for warm inflation, as will be seen in this section 
the simplest monomial models are consistent both from the cosmological and particle physics perspective. 

One generic feature of warm inflation models is that they require a large number of fields, usually in the hundreds 
or more. From the perspective of the simple single or few field inflation models typically seen in cold inflation, this 
feature of warm inflation models may appear undesirable. However in most high energy particle physics models, there 
are typically many fields present. In fact, from the perspective of string theory, where there are a huge number of 
fields, warm inflation models can look very compelling. Thus, in its own right, the requirement of a large number of 
fields distinguishes warm inflation models from their cold counterparts, but one can not say that they are necessarily 
more peculiar. 

To calculate a warm inflation model it requires solving the evolution equation (j2.13p and constraining the result 
with the density fluctuation amplitude and spectral index. In addition, for a first principles model, the dissipative 
coefficient T has been computed from an underlying quantum field theory model. One assumption in such calculations 
is that the microphysical dynamics determining T is operating at time scales much faster than that of the macroscopic 
motion of the inflaton background field and the expansion scale of the Universe. To realise such an adiabatic regime 
leads to a set of consistency conditions [Tsj . 



Tri>0/0,i7, (7.1) 

where here represents all relevant decay widths of the fields responsible for dissipation. 

Another challenge in realising warm inflation models from first principles is the effect of radiative and thermal 
corrections to the effective potential, which if too large would ruin the necessary flatness of the inflaton potential. 
There are competing requirements in that large dissipation prefers large couplings whereas controlling radiative and 
thermal corrections requires small couplings. Supersymmetry provides a means to achieve both these requirements. 
The observation is that SUSY will cancel local radiative corrections such as the loop corrections to the effective 
potential. Thus the use of SUSY models with interactions like the ones shown in Sec. |Vl However SUSY is ineffective 
in cancelling time nonlocal loop effects from which dissipation emerges. Of course there are still limitations. For 
one thing, since inflation requires a nonzero vacuum energy, SUSY must be broken, thus cancellations of radiative 
corrections are never perfect. Moreover, SUSY is also broken at finite temperature and so thermal loop corrections 
do not cancel exactly, although a significant amount of cancellation does occur [9^ . 

With these considerations in mind, several warm infiation models have by now been constructed. The first of these 
were the distributed mass (DM) models [lB|. In these models there are a set of bosonic fields Xi coupled to the inflaton 
field through shifted couplings. The interaction term in the Lagrangian which realizes such shifted couplings has the 
form, 

^l{<P-M,fxl (7.2) 

so that when (j) ~ Mi, the mass of the Xi field becomes small. In particular the mass is meant to get below the 
temperature scale in the Universe, so that these x fields can become thermally excited. Once thermally excited, as 
the background infiaton field evolves, it is able to dissipate energy into these fields. This results in a dissipative term 
in the infiaton evolution equation [l^ . If these mass scales are now distributed over a range of values that (f) will 
go through, then during evolution of 4>, some subset of these x fields will be light and generate a dissipative term. In 
order to control the radiative corrections in this model, SUSY has to be implemented. A simple superpotential that 
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realises this model is 



(7.3) 



Here the bosonic part of the chiral superfield $ = cj) + 9i/j + 6'^ F , with 9ip = d"ipa and 61 = O^Oa, is the inflaton field (j), 
and it interacts with both the bose and fermi fields of the chiral superfields Xi = Xi +^'0xi ~^^^^Xi ■ '^^^ potential terms 
of the Lagrangian are obtained from Eq. ()7.3p by standard procedures; the potential is Ly = J d'^xcP0W{^, {^i}) + 
h.c. and the auxiliary fields F and F^ are eliminated through the "field equations", dW/dF = dW/dF^. = 0, which 
results in the Lagrangian only in terms of the bose and fermi fields. This leads for the above superpotential Eq. (|7.3p 
to a A(/)^ inflaton potential with interactions to the Xi fields similar to Eq. (|7.2[) and corresponding interaction terms 
to the ip^. fermi fields. The mass scales are distributed along the line which traverses during the infiationary 
period. Moreover the (f>^ coupling must remain small for successful infiation, and in this SUSY theory it occurs since 
the renormalization group equations for the quartic coupling are proportional to the coupling itself, thus even if there 
is another large coupling this will not lead to a problem. This model can generate warm infiation with adequate 
e- foldings to solve the horizon and flatness problems 15] as well as produce observationally consistent primordial 
fluctuations . It should be noted that the most general superpotential would include a term in Eq. (|7.3p linear in 
the Xi fields, ^^Xi and this term has been eliminated by hand. This term induces a dependent mass term to all the 
Xi fields and so must be very small for the success of this model. The stability of the SUSY theory under radiative 
corrections allows this term to be eliminated by hand. However one can prohibit the linear term in the superpotential 
more elegantly by imposing a charge under some, for example GUT, symmetry so that these Xi fields are not singlets. 

In [93, [9^ the DM model has been shown as arising from a fine structure splitting of a single highly degenerate 
mass level. For typical cases studied in [HI, U^, it was shown in p^] that for significant expansion e-folding, > 60, 
if M « g\Mi+i — Mi\ denotes the characteristic splitting scale between adjacent levels, warm inflation occurred in 

the interval lO^M ^ ^ 3 x lO'^M and of note, at temperature M ^ T and not T at the much higher scale of 
the mass levels ^ IQ^M. What makes these massive states light is precisely the shifted mass couplings. In the 
string picture, this arrangement corresponds to a flne structure splitting of a highly degenerate state of very large 
mass, ~ M5, with the flne structure splitting scale several orders of magnitude less than the mass of the state, say 

M - Mgut ~ IQ-^Ms. 

In Ref. [98| the following string scenario was suggested. Initially in the high temperature region, some highly 
degenerate and very massive level assumes a shifted mass coupling to (j). Since all the states in this level are degenerate, 
at this point they all couple identically as g^ Tlii{4>~ ^Yxt- ^^'^ string then undergoes a series of symmetry breakings 
that split the degeneracy and arrange the states into a DM model ~ ^iYXd with < {Mi — Mi+i)/Mi <^ 1. 

This string scenario has several appealing features: 

(i) Strings have an ample supply of highly degenerate massive states. 

(ii) The generic circumstance is that as temperature decreases, many of the degeneracies will break at least a little, 
and for warm inflation a little is all that is needed. Moreover, warm inflation occurs when T is at or above 
the fine structure splitting scale but much below the scale of the string mass level. Thus, for the respective 
mass level, warm infiation is occurring in a low temperature region. This further supports the expectation that 
degeneracies for that level have broken. 

(iii) The shifted mass coupling to (j) is much more probable to occur to a single mass level, albeit highly degenerate, 
as opposed to the coincidence probability to several mass levels. 

(iv) There are minimal symmetry requirements for interactions. Since zero modes and any higher mass level modes 
fall into representations of the gauge and Lorentz groups, the interacting fields must tensor together to form 



The distributed mass models are the only warm infiation models constructed in which the fields directly interacting 
with the infiaton are thermally excited. In general it is too difficult to control adequately the thermal loop corrections 
to the infiaton effective potential to maintain the needed fiatness of the potential. This has led to the development of 
the two stage dissipative mechanism of warm infiation 16] in which the inflaton (/) is coupled to a set of heavy fields x 
and V'x which in turn are coupled to light fields y and ipy. The main point is the heavy fields arc not thermally excited 
and so the loop corrections to the inflaton potential are only from vacuum fluctuations, which can be controlled by 
SUSY. A generic superpotential that realises the two stage mechanism is 



singlets. 



Iccay 




(7.4) 



i=i j=i 
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where ^ = 4> + ipB + 9'^F, X = x + dip^ + d^F^ and Y = a + ei)„ + 6'^F„ are chiral superfields. The field will be 
identified as the inflaton in this model with (j) ~ (j}^ + t] and ((/>) — (f>c- In the context of the two stage mechanism 
X is the heavy fields to which the inflaton is directly coupled and these fields in term are coupled to light Y fields. 
In order to see the efi^ect of this interaction structure on radiative corrections to the inflaton potential, a particular 
model has to be chosen. Thus, considering the case of a monomial inflaton potential and adding the superpotential 
term W4, = \/A$3/3 so that 



W = W4, + Wi, 



this model generates at tree-level the inflaton potential 



(7.5) 



^o(0c 



(7.6) 



When ipc 7^ 0, there is a nonzero vacuum energy and so SUSY is broken, 
between the x V'x SUSY partners with in particular 



This manifests in the splitting of masses 



m 



^(3' - \'^9)4>l + V2mg(j), + Aw? 



'■■Ax 



The one loop zero temperature effective potential correction in this case is 

,2 \ 



1287r 



9 ^ 2 ,4 /, "^V-x 



2 « Vq( 



VXg(j)l , 
VXgcbl 



(7.7) 



which is further suppressed than the tree level potential Eq. (|7.6p and so will not alter the flatness of the inflaton 
potential. There are several first principles warm inflation models that implement the two-stage mechanism Eq. (|7.4p , 
which will be summarized here. 



A. monomial potential 



The general form of the inflaton monomial potential to be studied is, 



Vi^)^Vo{-^ , (7.9) 
\mpj 

with n > 0. Without enough dissipation, i.e., either for cold inflation with Q = 0, or only weak dissipation with 
Q < 1, where Q is defined in Eq. ()2.17p . these kind of models lead to infiation only for values of the inflaton field 
larger than the Planck mass mp. On the other hand, in the strong dissipative regime due to the larger friction term, 
slow-roll conditions Eqs. (|2.16p can be fulfilled for values of the field well below the Planck scale. Thus Eq. (|7.9p 
can be regarded from the effective field theory point of view, with the potential well defined below the cut-off scale 
nip; higher order term contributions suppressed by mp will be then negligible, without the need of fine-tuning the 
coefficients in front. In (4^ A/" = -^x-^Jecoy ^^'^ 9* "were treated as free parameters, and it was examined for which 
values r]/Q and T/(j) can be kept small enough for at least 50 e-folds or so (and T/H > 1). For example, for a quartic 
potential with T^(0)^/^ ~ O.Smp and (f){0) = mp, in order to satisfy all the constraints it required 5* < 100 but 
Af > 2300. Similar results were obtained for other powers of the potential. By lowering the value of the potential, it 
was found to be easier to fulfill all conditions except that for the ratio rj/Q. Keeping the latter below one gives the 
lower bound: 

.3/4_„ , . .7/4 



> 8.4 X 10-2 Li/7(„ _ 1) _^ ^^1 ^ (7 



T/(0)i/4 



7 



and the lower V^(O) is, the larger has to be. For example, for n = 4, V{Oy^^ /mp = 0.1, and = 10, it requires 
Af > 2800, but getting to the number of degrees of freedom for the MSSM, 5, = 228.75, would require A/" > 29000. 
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The interesting result from this analysis was that due to the extra friction, inflation occurred for values of the field 
below the Planck scale, although the model prefers an initial value of the height of the potential only an order of 
magnitude or so below the Planck scale. 

On the other hand, the amplitude of the primordial spectrum is also affected by the strong dissipative friction 
term and the presence of a thermal bath. In order to keep the amplitude of the primordial spectrum consistent with 



WMAP's value M 



)l/2 



5.5 X 10 ^, it required a potential much smaller than 0(10 



3). For such a value of 



the potential, it needs roughly Af ^ O(IQ^) in order to get at least 50 e-folds in the strong dissipative regime. 



B. hybrid potential 



In 143 small field models of inflation were also considered of the form, 



VicP) - Vo ( l + /nn^ 



71 > 



71 = 



(7.11) 



(7.12) 



Given that during inflation the potential is dominated by the constant term Vq , the value of the field can easily be kept 
below the Planck scale in th ese models during slow-roll inflation. These potentials can be regarded as a generalisation 
of a hybrid model [isl. Il00j |. where inflation ends once the inflaton field reaches the critical value, destabilising the 
waterfall field coupled to it. Those interactions are not relevant to study the slow-roll dynamics, only to mark the end 
of inflation, and therefore they do not need to be considered in the inflationary potential Eqs. (|7.1ip and (|7.12p . As 
observed in ^431 , the same interactions between the inflaton and the waterfall field required by the hybrid mechanism 
will give rise to dissipation, and leads in the low-T regime to the dissipative coefficients given in Subsec. IV Bl The 
case 71 = 2 would be the standard hybrid model [H, llOOj ] , with a mass term for the inf laton, wh ereas 71 = is the 
susy model with the logarithmic correction coming from the one-loop effective potential |lOll . llO^ . 

In supersymmetric hybrid models, one needs to worry about the 77-problem discussed earlier |48l. |49| . [sol . [sil . [s^ . 
i.e., the fact that generically SUGRA corrections give rise to scalar masses of the order of the Hubble parameter, 
including that of the inflaton, which in turn forbids slow-roll inflation. Different solution to this probl em exist i n th e 
literature, for example by combining specific forms of the superpotential and the Kahler potential [48l . ll0lLll02l . llOSt . 
Nevertheless, typically, although the quadratic correction can be avoided, i.e., a mass contribution, SUGRA corrections 
manifest as higher powers in the inflaton field [47] . In the case of strong dissipative warm inflation, the presence of 
the extra friction term alleviates the problem: slow-roll conditions are fulfilled also for inflaton masses in the range 
H ~ 7710 ^ V HT. In addition, the values of the field being smaller than in standard cold inflation, the effect of higher 
order SUGRA corrections is also suppressed. 

It also should be noted that in all the warm inflation models constructed in this Section, the inflaton field <I> is 
a singlet, for which an extra complication arises in that a linear term in $ is allowed in the Kahler potential. This 
induces a tadpole term for the singlet, which results in a large vacuum expectation value (VEV) of the order the cutoff 
scale. This can lead to problems for both the low-energy theory in destabilizing the electroweak scale jl04l | as well as 
make it difficult to realize inflation [10511 . Slight modifications have been shown can keep the singlet corrections under 
control for both the low-energy theory [106, 107] as well as inflation jl05l . ll08| . Without a higher theory, the terms in 
the Kahler potential ultimately are arbitrary and one can always tune the linear term to be small. Alternatively one 
can impose a symmetry on $ which prohibits such a term. As yet, warm inflation model building has not explored 
these detailed questions. However one point should be noted, that the presence of the friction term in warm inflation 
allows for larger 77 and e, and so the effect of all terms including the linear term in the Kahler potential is alleviated. 

The amplitude of the primordial spectrum is given by Eq. (2.28) and the spectral index 715 is found in [43| to be 



and 
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(7.13) 



(7.14) 



where it is assumed e <^ rj. Notice that the spectral index is of the order of 0{ri/Q), whilst the running is of the 
order of 0{rf' /Q^). Therefore, the same condition needed to have slow-roll in the strong dissipative regime will avoid 
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having too large a spectral index. The model has a blue-tilted spectrum when n < 7, including the case of n = 0, i.e., 
the logarithmic potential. The more negative running, the more blue-tilted the spectrum can be, which would be the 
case for < n < 7 with n'g « —{ns — l)^n/(21 — 3n). 

Given that the field decreases during inflation, so does rj/Q, and also pr/V (or equivalently T/H) for any power 
n ^ 0, being constant for the logarithmic potential n — 0. Therefore, the energy density in radiation will never 
dominate in this regime. On the other hand, (j)/T diminishes for n < 4, but Q only decreases for n > 2. Therefore, 
for a logarithmic or quadratic potential once the system is brought into the strong dissipative regime, it stays there 
until the end. Indeed one can start in the weak dissipative regime, with T > H but Q < 1, and it will evolve into 
Q > 1. 

For n = the model exhibits the interesting situation of a transition from "cold" —^ "weak warm" — > "strong warm" 
inflation. For n ~ 2, T remains constant until it diminishes in the strong dissipative regime, so that the parameter 
space divides into either cold or warm inflation, but the system can evolve from weak to strong dissipation. In the 
weak dissipative regime, quantum fluctuations of the inflaton field will also have a thermal origin, with spectral index 

[loi] 

ns ~l-2e + 2^^, (7.15) 

n — 1 

where note that when n = 1, it also follows rj = and so there is no singularity. Notice that even in the weak dissipative 
regime a logarithmic potential gives rise to blue-tilted spectrum, while n > leads to a red-tilted spectrum, just the 
reverse than the standard cold predictions. For larger powers n > 2, the opposite behavior is found, so even if inflation 
starts in the strong dissipative regime it will evolve towards the weak and the cold regime. When this happens before 
the last 50 e- folds of inflation, then dissipation becomes irrelevant. 

There are a few specific interesting features of these models found in [42], which are discussed next for n — 0,2, 4. 

Case n = 0: Hybrid logarithmic potential. This model has its lower bound on JV in the limiting case for 
slow-roll warm inflation, when T/H ~ 1, Q{0) ~ 1 and {(f) /T)n — 10, which gives: 

A/" - 11.87 exp[-0.23iVe.gyVA/'^^^] • (7-16) 

For example, with .g* ~ 10, = 50 we have the lower bound: TV ~ 180, (/)(0)/i/ ~ 150, ~ 244, 'q/Q ~ -0.054; 
with ~ 228.75, Ne = 50 we have: M ~ 1350, 0(O)/iJ ~ 1.2 x 10^, ~ 5.2 x 10^ and 77/Q ~ -0.1. Those values 
oi N = N^N^ecay quitc in the range of a realistic model; for example with x in the 126 or 351 of S'O(IO) {Eq), 
and N'^i^^^y ~ C(10), one can expect having M in the range of a few thousands. However, for such values in [i^, the 
amplitude was found to be too large. In order to match the amplitude of the primordial spectrum with WMAP's 
value, it requires a larger initial value of the field (/)(0), but then the value of {M) has to be larger in order to stay 
within the low T approximation, with (p/T > 10. Satisfying WMAP's constraints requires then Af > 5 x 10^, which 
is rather large, with 0(O)/i? > 2 x 10^ and a-^ > 8 x 10^^, where 

'nVo\ fH' 



0. 



An alternative example in p3 | was inflation in the weak warm regime, and a transition from weak to strong 
dissipation at the end. This helps in fixing the amplitude of the spectrum to lower values. The following conditions were 
still imposed, that (a) T/H > 1 (b) to obtain enough inflation, i.e., Ne « 50. These translated into Af > O.OSt/^/ry^, 
with r] < l/{2Ne), and so gives the lower bound TV > Q.2g^N'^. For example for = 288.75 and — 50 it gives 
TV > 1.2 X 10^, which again is rather large. 

Case n = 2: Hybrid Quadratic potential. Figure ^ shows results for the evolution of the ratios Q, T/H and 
T/(j), for TV = A^x^Jecay = 10000, 20000, 30000 and = 228.75. The end of inflation was taken as the point at which 
T = H, and the analytical approximations do not any longer hold, and therefore in the figure counts the number 
of e-folds left to the end of inflation. The condition T > H \n turn translates into a lower bound for TV: 

N > 0.052g, + ' . (7.18) 

From the approximated expression for the spectral index, Eq. (|7.13p . if one wants to keep ns within the observable 
range, it requires ?//(3o < 0.093, which for ~ 50 gives the lower bound TV > 32.5(7*. Again, having slow-roll 
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FIG. 6: Warm inflation for hybrid quadratic potential: Evolution of the ratios Q (top lines), T/H (middle lines) and T/(j) 
(bottom lines) depending on the number of e-folds to the end of inflation, for different values of A/'=10000, 20000, 30000, with 
= 228.75, (f>{Q)/mp = 0.21, 77 = 3, and V^'^ /mp = 3 x lO"*. 

warm inflation for example with ~ 228.75 needs M > 7500, but for 3* ~ 10 it only requires M > 325. As an 
example, Fig. [7] shows the predicted spectral index depending on the number of e-folds left to the end of inflation, for 
Af = N^Nj^^^y = 10000, 20000, 30000 and 5* = 228.75. The corresponding spectral index of the primordial spectrum 
would be that at around 50-55 e-folds, which is always ns < 1.2. The value of the running can be obtained from Eq. 
(fm)) . and it is given respectively by n'g ~ -2.5 x 10'^, -8.5 x lO""*, -4.7 x 10"''. 

Controlling the amplitude for the primordial spectrum from getting too large was found to require values of 77/Qo 
as large as possible, but not too large values of Qo- For values of A/", 5* within the range of Eq. (|7.18p . the amplitude 
remains below say 10"'' for values of Qo of order 0(10). Therefore in these kinds of models parameter values can be 
found giving rise to the right order of magnitude for the primordial spectrum in the strong dissipative regime, but 
the stronger constraint comes from avoiding an overly blue-tilted spectrum. The results for this example also reveal 
one generic feature of warm inflation solutions in quantum field theory models, that the total duration of inflation 
tends to be small, of order the observational requirement of « 60 e-folds. This is in contrast to cold inflation models, 
which in general cases can predict huge numbers of e-folds, orders of magnitude larger that the observational lower 
bound. This fact about smal l e-fo lds in warm inflation has been explored as a possible solution to the low quadrapole 
observed in the CMB data Moreover, the small e-folds predicted in warm inf lation cou ld also be a possible 

benefit to the transplankian problem, in which small e-folds of inflation are preferred [llll Ill2j | . 

Case n > 4: Hybrid quartic and higher powers. In this case having 50 e-folds in the strong regime requires 
for example N = N^N^ecay ^ for g.^. ~ 10, and TV > lO'^ for ~ 228.75. In addition if we want to get the 
right amplitude for the spectrum and spectral index, it increases Af by one order of magnitude, since the quartic 
coupling ~ \ has to be adjusted to rather small values, and then the initial value of the field to larger values to 
have 0(O)/r > 10. Numbers do not change much whether 10 or 50 e-folds of inflation are demanded in the strong 
dissipative regime. 
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FIG. 7: Hybrid quadratic potential: Spectral index depending on the no. of e-folds to the end of inflation, for different values 
of 7V=10000,20000,30000, with 5, = 228.75, (fi{0)/mp = 0.21, = 3, and Vo^^/mp = 3 x lO""*. 
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FIG. 8: The figure shows the log C^-log |r;| plane, where = 0.64:h'''NxNl^^^y. The light grey (light green) area represents 
the space where the thermal spectrum of perturbations matches the observed amplitude A^* ~ 50 e-foldings before the end of 
inflation, with ^ond < 10 mp. Within the dark grey (reddish) band, the resulting value of the spectral index is within the 1-a 
window; n = 0.960j;o;oi3, as inferred by WMAP+BAO+SN data in the ACDM model for negligible tensor perturbations. In 
the graph lines of constant H and Q = T /3H axe also depicted. 



C. hilltop potential 

In [3| the warm hilltop model [9lj was investigated for the potential 

V = Vo-^\m^\^^ + ..., 

where Vq = iH^m?p, m? — V"(0), and the dots represent higher order terms that become important only after relevant 
scales exit the horizon during inflation. The inflaton fleld <f) is coupled to the flelds of the two stage mechanism in 
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Eq. (|7.4p . The scenario starts with the inflaton field close to the hilltop. In this model was constrained to obtain 
adequate e-folds of warm inflation and a consistent amplitude for density perturbations. In addition constraints were 
placed to avoid gravitino overproduction. The resulting parameter space for the strong dissipative warm inflation 
regime is shown in Fig. [51 where = O-dAh'^N-^N^^^^^. Moreover nongaussian effects in warm inflation [56llll3llll4l | 
hav e bee n studied for this model in [4^ . In the strong dissipative regime, there are in general large nongaussian effects 
[s^ . Ill4| . In particular it was shown in [sBl that entropy fluctuations during warm inflation play an important role in 
generating non-Gaussianity, with the prediction 

-151n(l + |)-5</„.<f,„(l + i|)-^, (7^9, 

where /nl is the non-linearity parameter and r = T/5H. For the warm inflation results in Fig. [51 Q ranges from 
10 to 10^, and this implies from the above equation that |/nl| ra nges from 10 to 180. This is an interesting result 



in light of the recent WMAP analysis: the third year CMB data |115{ gives 26.9 < /nl < 146.7 at 95% confidence 
level, although the five-year WMAP data ^ give the limit -9 < /nl < HI (95% CL). The latest data show then 
a tendency for /^v^ > 0, alth ough this still will need to be confirmed by future data, and in particular by data 
from Planck surveyor satellite |ll6l |. If this is the case, this would disfavor conventional cold inflation models which 
generally yield very low values of /nl ^ 1- On the other hand, the strong dissipative warm inflation regime, such as 
the one found in this hilltop model, would be consistent with a non-gaussian signal. For example, in this case taking 
/nl ^ 110 translates in an upper limit on Q < 1.3 x 10^, and therefore from Fig. [5] on a lower limit on the scale of 
inflation given hy H > 10® GeV. 

VIII. CONCLUSIONS AND FUTURE PERSPECTIVES 

Generically, the inflaton interacts with other flelds in any typical inflation model and so its dynamics is dissipative. 
As such, inflation, like most dynamics in nature is an open system phenomenon, thus requiring a much more complex 
analysis of its dynamics than the one typically formulated for the cold inflation picture. This general point has 
been voiced by B. L. Hu and coworkers [71| and more specifically in the initial motivating papers of warm inflation 
[1, [l^, [l3|- In all these cases, the point has been made that the problem of inflation encompasses many different 
branches of Physics, from nonequilibrium statistical dynamics to particle physics phenomenology. 

Warm inflation dynamics is a rich area of study for both quantum fleld theory real time dynamics and for particle 
physics model building. The study of the warm inflation dynamics has almost exclusively motivated the understanding 
of strong dissipative behavior in quantum field theory. This started initially with the work of Berera, Gleiser and 
Ramos [13|, in which dissipation was examined in the overdamped regime for the first time in quantum field theory 
using extended linear response calculations. Since then, specific progress has been made to underpin interaction 
structures in interacti ng q uantum field theory models which lead to strong dissipative behavior under warm inflationary 
conditions [H, [t^, [H, iSSll. A t a more general level, this work has motivated the first calculations of dissipation in the 
low temperature regime In this review, a physical picture has also been developed in Sec.[Vl for explaining the 
dissipation behavior found in warm inflation from the quantum field theory calculations, and this is further developed 
in [84]. 

Up to now all these studies has been based on variants of linear response methods, in which case the dynamics 
is supposed to happen close to equilibrium, or in a quasi-adiabatic regime, including various types of resummations. 
Several extensions to this work, which will give a more accurate understanding of dissipation, are under way. For 
instance, when departing from the quasi-adiabatic regime for the field dynamics, it is expected that the local Markovian 
approximation commonly used to analyse the field equations can differ significantly from the exact nonlocal equations. 
Preliminary tests have shown that this difference can be very large for a period of time starting from the initial period, 
but with differences between the dynamics getting smaller at longer times, depending on the model parameters [117|. 

The methods used in this review are all based on the effective equation of motion derived from the action functional. 
In order to study strong nonequilibrium dynamics, it requires methods beyond these. In this case, full kinetic set 
of equations for the relevant fields have to be studied, which contain information not only of the inflaton effective 
dynamics but also about thcrmalisation and equilibration. Work in this direction, also can help to better understand 
the physics of particle production during the system dynamics and be used to check the reliability of using the 
approximation of thermal initial conditions for the bath flelds. Initial work in this direction, though not dir ectly 
related to the type of models relevant to warm inflation as discussed extensively e.g. in Sec. IVIIl has appeared [8l| . 
while a work more on the warm inflation dynamics motivated side is under way [ll8|. 

This review has presented in detail the methods used so far to understand warm inflation dynamics as well as the 
limitations of these studies. In particular these methods used so far are primarily quasi-adiabatic approximations 
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for the fields with the assumption of near thermal equilibrium evolution. Despite these limitations, these results find 
parameter regimes in which physically acceptable solutions exist over time periods sufficiently long to be of use in 
studies of warm infiation. During this time interval in which these approximations apply, and where a local Markovian 
dynamics can be used, in contrast to the full non Markovian one, in the typical model implementations discussed 
in Sec. IVIIi the amount of radiation production was seen to be sufficient to change the cold inflationary picture 
predictions regarding the density perturbations, thus requiring its description in terms of the warm inflation picture. 
Moreover, once dissipative effects are strong enough, inflation can be sustained and driven longer than when these 
effects are neglected. Consequently, parameter values typically required in the cold inflation case can be relaxed, 
which can help to evade various problems that plaque the standard scenario of inflation, like the graceful exit and 
77-problem, discussed in details in this revie w, as well as the problems of quantum-to-classical transition [ol. Ill9l | and 
the initial conditions for inflation [l20l . Il2l( . 

The development of the quantum field theory dynamics of warm inflation has in turn been applied to particle physics 
model building, in which warm inflation dynamics in realised. Early on it was recognisied in [ll| that warm infiation 
has some appealing and unique model building features. In particular, it offers a simple solution to the 77-problem and 
for monomial potentials, observationally consistent inflation can occur for the infiaton amplitude below the Planck 
scale (0) < mp. These features have been realised in explicit flrst principles quantum field theory models of warm 
infiation in [ll|, [l^ [Z^, . These studies are now being extended to develop a complete particle cosmology in which 
not only is warm infiation realised, but in addition other features such as leptogencsis and gravitino abundances are 
addressed, developing in depth some of the work already started on these topics 42, 44, 109, 122, 123]. In a sep arate 
direction, the warm infiation models developed so far in [i^, |3| have been in the low temperature regime. In jl24| 
this is being extended to higher temperatures. This requires calculating all the thermal loop corrections in the SUSY 
models that have the two-stage interaction structure Eq. (17. 4p relev ant for warm infiation. Some initial work has 
been done in , and a more detailed analysis is now underway [l24l | . 
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